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Хімічна  взаємодія  атмосферних  викидів  речовини  з  компонентами 
атмоефери  породжує  нові  еполуки  й  речовини,  дія  яких  на  живі 
організми,  рослинний  світ,  метали,  споруди,  будівлі  може  мати 
деетабілізуючий  вплив,  бути  шкідливою,  викликати  у  організмі  людини 
незворотні  біологічні  процеси.  Боротьба  з  викидами  шкідливих  і 
токсичних  сумішей,  екепериментальні  доелідження,  фізико-математичне 
моделювання  та  розробка  фізико-хімічних  способів  конверсії  сумішей  є 
наукоємними  етапами  і  потребують  викориетання  значного 
інтелектуального  рееуреу. 

До  найбільше  ефективних  методів,  що  вирішують  проблеми 
забруднення  навколишнього  середовища,  відноеяться  електрофізичні  - 
заетоеування  електричних  розрядів  та  електронних  пучків,  іонізуюче 
випромінювання,  плазмоетруйні  реактори  [1-8].  До  переліку  цитованих 
робіт  слід  додати  нашу  роботу,  у  якій  вперше  була  розв’язана  задача 
динаміки  елементарного  акту  хімічної  реакції  на  твердій  поверхні  [9].  У 
цій  роботі  на  прикладі  молекули  СО  показаний  перехід  хімічного 
зв’язку  від  етану  макеимальної  етійкоеті  до  моменту  його  дисоціації  у 
процесі  наближення  до  точкового  електричного  заряду  (кулонівського 
центру),  у  дослідженнях  впливу  електричного  точкового  заряду  на 
етійкість  хімічного  зв’язку  вивчаєтьея  як  мінімум  два  питання:  перше  - 
це  залежніеть  руйнування  хімічних  зв’язків  від  щільності  розподілу 
зарядів  на  поверхні  або  концентрації  у  проеторі  газового  еередовища, 
відетані  молекули  від  заряду  та  температури;  друге  -  утворення  нових 
зв’язків  та  нових  речовин  в  залежності  від  термодинамічних,  кінетичних 
та  інших  факторів. 

В  роботах  [10-13]  доеліджувалися  процееи  дисоціації  водню  у 
наносекундному  розряді,  релаксаційні  та  плазмохімічні  процеси  після 
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наносекундних  газових  розрядів.  Під  час  електрофізичних  обробок 
речовини  виявлено  ряд  нових  закономірноетей,  які  слід  ураховувати  при 
переході  від  екепериментів  до  практичного  викориетання.  Ветановлено, 
наприклад,  що  зменшення  тривалості  розряду  збільшує  за 
експоненціальною  залежністю  ефективність  синтезу  нітридів  фосфору  і 
конденсованого  озону.  Детально  досліджена  конверсія  оксидів  азоту  і 
сірки  в  продуктах  згоряння,  ініційована  стримерною  короною, 
вивчаються  процеси  коронування  катодів  в  електричних  розрядах, 
досліджуються  процеси  деконтамінації  (знезаражувальня)  повітря  і  таке 
інше. 

Одним  із  напрямків  вирішення  проблеми  очистки  повітря  від 
домішок  сірковуглецю  є  створення  в  газі  слабкоіонізованої  плазми 
електричними  розрядами  або  електронними  пучками  [1;2;8].  За  рахунок 
підбору  параметрів  пучка  і  заряду  одержують  оптимальний  режим 
конверсії  токсичної  домішки.  Експериментально  встановлено  [1]  два 
режими  конверсії  домішки,  які  залежать  від  типа  використовуваного 
пучка,  відрізняються  механізмами  та  кінцевими  продуктами:  у  першому 
випадку  продуктами  конверсії  являються  оксиди  елементів  газової 
суміші,  у  другому  -  сполука  (-С8-)п  у  конденсованому  стані.  Слід 
відмітити,  що  у  плазмі,  створеній  пучком  електронів  високої  щільності 
дія  зовнішнього  електричного  поля  ініціює  утворення  оксидів  азоту  і 
зниження  інтенсивності  конверсії  сірковуглецю.  Інший  еценарій 
складається,  якщо  впливати  на  суміш  слабкострумовими  пучками.  При 
збільшенні  вкладу  іонно-молекулярних  процесів  (головним  чином  за 
рахунок  зниження  інтенсивності  електрон-іонної  рекомбінації) 
формуються  кластери  і  у  них  відбувається  полімеризація  сірковуглецю. 
Тут  слід  підкреслити,  що  ефективність  видалення  молекул  сірковуглецю 
із  іонізованого  повітря  стає  дуже  чутливою  до  слабкого  електричного 
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поля.  Як  вважають  автори  [1]  це  пов’язано  з  залежністю  швидкості 
генерації  іонів,  утворюючих  кластери,  від  напруженості  поля. 

Особливої  актуальності  набувають  питання  безпеки  середовища  та 
надійності  роботи  електроустаткування  у  підземних  виробках. 
Вибухозахищене  і  рудникове  електроустаткування  повинне  забезпечити 
безпеку  застосування  електроенергії  і  високу  надійність  у  специфічних 
умовах  підземних  підприємств  гірничодобувних  галузей  промисловості 
(шахт,  рудників),  а  також  підприємств  хімічної,  нафтової,  газової  та 
інших  галузей  народного  господарства,  що  мають  вибухонебезпечну 
атмосферу.  У  той  же  час  підземні  електроустановки  експлуатуються  в 
умовах  підвищеної  вологості,  запиленості  і  вибухонебезпечності 
рудникової  атмосфери. 

Фізико-хімічні  процеси,  виявлені  під  час  виконання  даної  роботи, 
з'явилися  фактично  відкриттям  для  фахівців,  що  розробляють 
вибухозахищене  електроустаткування,  і  визначили  головний  напрям  в 
створенні  такого  електроустаткування:  розробку  і  вживання 
комутаційних  апаратів,  які  здійснюють  гасіння  електричної  дуги  в 
герметизованих  камерах.  За  ініціативою  УкрНДІВЕ  (м.  Донецьк)  вперше 
в  Україні  було  організовано  виробництво  контакторів  з  вакуумними 
дугогасильними  камерами  на  струми  160,  250,  400  А  та  напруги  до 
1200  В,  на  базі  яких  випускаються  всі  сучасні  пускачі  та  станції 
управління. 


6 


1.  УМОВИ  ТА  ПАРАМЕТРИ  УТВОРЕННЯ  АГРЕСИВНИХ, 
ПОЖЕЖО-  ТА  ВИБУХОНЕБЕЗПЕЧНИХ  СЕРЕДОВИЩ 
У  ПІДЗЕМНИХ  ВИРОБКАХ  ШАХТ 


Для  збільшення  потужності  і  продуктивності  гірничих  машин  і 
механізмів  необхідно  підвищувати  напругу  живлення  цих  машин. 
Підприємства  вугільної  і  деяких  інших  галузей  гірничодобувної 
промисловості  перейшли  на  живлення  технологічного  устаткування 
напругою  660  Вив  теперішній  час  проводяться  роботи  з  масового 
переводу  систем  електропостачання  на  наступну  ступінь  напруги  - 
1140В.  Однак  досвід  експлуатації  вибухобезпечних  комутаційних 
апаратів  показав,  що  спроби  швидкого  переведення  підземних 
електроустановок  на  підвищену  напругу,  окрім  недосконалості  систем 
захисту,  були  затримані  на  декілька  років  ще  одним,  раніше  невідомим 
чинником:  катастрофічно  швидким  зниженням  опору  ізоляції  і  виходом 
з  ладу  вибухозахищених  і  рудникових  комплексних  апаратів  управління 
і  захисту,  які  працюють  в  інтенсивних  режимах  комутації  силових 
електричних  мереж  при  напрузі  660  В  і  вище.  Наприклад,  термін  служби 
пускачів  на  660  В  складав  від  декілька  тижнів  до  3-х  місяців. 

Атмосфера  підземних  гірничих  виробок  (температура,  вологість, 
загазованість  і  запилення)  впливає  на  стан,  термін  служби  і  надійність 
вибухобезпечного  електроустаткування.  Але,  додатково  до  зовнішніх 
впливів  гірничої  атмосфери,  в  обмеженому  обсязі  рудникових 
вибухобезпечних  апаратів  формується  особливий  мікроклімат, 
відмінний  по  своїх  параметрах  від  клімату  гірничих  виробок. 
Вибухонепроникна  оболонка  апарата  перешкоджає  вільному  обміну 
внутрішньої  мікроатмосфери  з  атмосферою  навколишнього  середовища. 
Окисли  азоту  з'єднуються  з  парами  води,  що  є  присутніми  у  вологому 
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повітрі,  і  утворюють  азотну  й  азотисту  кислоти.  Ці  активні  окислювачі 
осідають  на  контактні  з'єднання,  ізоляційні  матеріали,  металеві  вузли  та 
деталі,  що  призводить  до  окиелювання  контактних  і  рухливих  з'єднань, 
зниження  опору  ізоляції,  причинами  цього  є,  наеамперед,  відкрита 
електрична  дуга,  що  виникає  при  комутації  підключеного  до  апарата 
навантаження,  а  також  теплові  та  фізичні  процеси  усередині 
вибухонепроникної  оболонки  при  проходженні  етруму  навантаження  по 
еиловим  і  контрольним  апаратам  і  етрумопроводам. 

Для  створення  і  освоєння  виробництва  на  підприємствах  України 
нового  покоління  вибухозахищеного  електроустаткування,  що 
забезпечує  безпечне  і  безперебійне  електропоетачання 
високопродуктивної  гірничодобувної  техніки  для  підвищення 
продуктивноеті  праці  і  об'ємів  видобутку  кориених  копалин,  необхідно 
було  вирішити  цілий  ряд  доелідницьких,  конетруктореьких  і 
технологічних  задач,  одна  з  яких  полягала  у  ветановленні  і  доелідження 
причин  прискореного  руйнування  ізоляційних  матеріалів  і  покриттів 
металевих  деталей  вибухозахищених  апаратів  управління  при  напрузі 
660  В  і  вище,  а  також  в  розробці  методів  і  заеобів  усунення  вказаних 
негативних  явищ. 

В  результаті  доеліджень,  еуть  яких  представлена  в  даній 
монографії,  розроблена  фізико-математична  модель  розрахунку 
електронних  молекулярних  термів  з  врахуванням  екранування  ядерних 
потенціалів  електронними  остовами  зв'язаних  електронів.  Особлива 
увага  відводилася  розробці  методики  розрахунків  рівнів  коливально- 
збуджених  етанів  молекул  газу.  Відкрита  нова  закономірніеть  розвитку 
високошвидкіених  физико-хімічних  реакцій  в  рудниковій  атмосфері, 
обумовлена  об'ємними  і  поверхневими  процесами  в  нерівноважній 
низькотемпературній  плазмі  газового  розряду,  які  еупроводжуютьея 
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руйнуванням  молекул  та  утворенням  нових  з'єднань  під  впливом 
електричного  поля  точкових  зарядів. 

Високоенергетичний  вплив  є  одним  з  найбільш  ефективних 
способів  зміни  структури  речовини,  складу,  стану  і  багатьох  її 
властивостей.  Одним  з  таких  способів,  що  роблять  помітними  впливи  на 
стійкість  молекул  та  фізико-хімічні  властивості  повітря  у  підземних 
виробках  щахт  і  рудників,  є  нерівноважний  електричний  розряд  у  газах 
[14-18].  У  газовому  розряді  енергія  електричного  поля  передається 
електронам  слабо  іонізованої  плазми,  які  в  свою  чергу  забезпечують 
надтеплове  заселення  збуджених  станів  атомів  і  молекул.  Нейтральні 
молекули,  збуджені  електронами,  можуть  стимулювати  перетворення,  що 
здійснюються  у  нерівноважних  плазмохімічних  системах. 
Нерівноважність  може  бути  обумовленою  різними  фізичними  впливами  - 
електромагнітного  поля,  надзвукового  витікання,  а  також  самою 
хімічною  реакцією,  тому  що  хімічна  реакція  як  пороговий  процес 
зменшує  кількість  молекул  з  енергією,  що  перевищує  порогову,  і  змінює 
тим  самим  функції  розподілу  молекул  за  енергіями. 

Великою  різноманітністю  визначаються  механізми,  за  якими 
здійснюються  у  нерівноважних  умовах  хімічні  реакції.  Один  і  той 
сумарний  процес  у  залежності  від  ступеню  іонізації,  електронної 
температури,  тиску  та  питомої  енергонапруженості  розряду  може 
здійснюватися  принципово  різними  шляхами,  з  суттєво  відмінною 
енергетичною  ефективністю.  Механізми  плазмохімічних  реакцій  мають 
ряд  особливостей,  обумовлених  наступними  факторами: 

•  реакції  дисоціації  хімічних  зв’язків,  в  результаті  яких 
утворюються  вільні  радикали  (ініціюються  такі  реакції  коливально-  та 
електронно-збудженими  молекулами  і  електронами); 

•  дисоціація  хімічних  зв’язків  через  електронно -збуджений  стан, 
що  протікає  у  дві  стадії:  1  -  електронне  збудження,  2  -  дисоціація 
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збуджених  станів; 

•  іонно-молекулярні  реакції  за  участю  електронно-збуджених 

іонів. 

У  нерівноважних  системах,  в  тому  числі  електричному  розряді  в 
газах  кінетика  хімічних  процесів  суттєво  відрізняється  від  звичайної 
хімічної  кінетики.  Нерівноважна  хімічна  кінетика  враховує  квантову 
енергетичну  структуру  молекул  і  атомів,  тобто  концентрацію  кожного 
компонента  у  кожному  енергетичному  стані  та  канали  хімічних  реакцій. 
Взагалі  хімічні  реакції  у  низькотемпературній  плазмі  є  багатоканальними 
процесами  і  визначають  таким  чином  широкі  можливості  управління  їх 
напрямками. 

Важливу  роль  коливальне  збудження  молекул  відіграє  в  організації 
нерівноважних  процесів  і  селективних  реакцій.  Для  нерівноважної 
плазмохімії  механізм  коливального  збудження  молекул  електронним 
ударом  характеризується  спецефічною  перевагою,  яка  полягає  в  тому,  що 
цей  механізм  може  стимулювати  хімічні  перетворення  з  найбільшою 
енергетичною  ефективністю.  Іншими  словами  коливальна  ступінь 
свободи  молекул  є  слідством  дії  розрядного  енерговкладу.  Для  цілого 
ряду  неелектронегативних  молекул  газу  -  N2.  СО,  СО2,  Н2  та  інших  - 
коливальна  енергія  реагентів  є  найбільш  ефективною  у  подоланні 
активаційних  бар’єрів  елементарних  реакцій.  Відомо,  що  при 
коливальному  збудженні  молекул  енергетичний  поріг,  як  правило, 
співпадає  з  активаційним  бар’єром  реакції,  але  через  електронне 
збудження  молекул  виявляється  помітно  нижчим  за  поріг  самої  реакції. 
Наприклад,  коливальне  збудження  молекули  водню  потребує  для 
дисоціації  4,4  еВ  (енергія  зв’язку),  а  дисоціація  електронним  ударом 
потребує  у  двічі  більшої  енергії.  Для  інших  молекул  приклади  аналогічні. 

Як  один  із  прикладів  виникнення  хімічно  активної  плазми  з 
нерівноважним  коливальним  збудженням  молекул  газу  спостерігається  у 
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просторі  МІЖ  електродами  під  чає  розімкнення  контактів,  що  нееуть 
електричний  струм,  та  на  поверхні  аноду  як  ефективного  джерела 
вторинних  атомів,  що  утворюютьея  під  дією  потоку  швидких  іонів  [15- 
18].  Утворення  вторинної  плазми  обумовлено  проникненням  в 
міжелектродний  проміжок  вторинних  атомів,  що  іонізуються. 
Дослідженнями  ветановлено,  що  збудження  і  дисоціація  молекул  газу 
здійенюєтьея  не  тільки  під  час  електричного  розряду  безпосередньо,  але 
й  піеля  іенування  плазми. 

Між  електродами  у  газі  акти  дисоціації  хімічних  зв’язків  та  хімічні 
реакції,  що  протікають,  можуть  бути  ініційовані  головним  чином  за 
такими  оеновними  умовами  та  причинами:  в  плазмі  під  чає  та  унаелідок 
газового  розряду  енергія  електричного  поля  передається  електронам 
слабо  іонізованої  плазми,  які  в  свою  чергу  забезпечують  надтеплове 
заеелення  збуджених  етанів  атомів  і  молекул;  в  плазмі  на  поверхнях 
катоду  та  аноду;  за  рахунок  впливу  на  хімічний  зв’язок  напруженоеті 
поля  точкових  зарядів,  розташованих  на  поверхнях  [9;  15;  16]. 

Результати  доеліджень  в  облаеті  фізики  і  хімії  електричних 
розрядів  у  газі  добре  відомі  та  мають  іеторію,  що  екладаєтьея  майже  із 
двох  етоліть.  До  фундаментальних  результатів  слід  віднести  механізми 
утворення  і  руйнування  збуджених  атомів,  явища  переноеу  та  кінетика 
збуджених  молекул  у  плазмі,  закономірноеті  елементарних  радіаційних 
процееів,  переноеу  випромінювання  у  еередовищі,  механізм  переноеу 
електричного  заряду  і  багато  інших  фізичних  ефектів.  Що  стоеується 
наукових  доеліджень  в  облаеті  взаємодії  хімічних  зв’язків  молекул  з 
точковими  зарядами,  то  таких  робіт  дуже  мало  і  вони  обмежуютьея 
декількома  десятками  наукових  статей. 

Утворення  холодної  нерівноважної  плазми  при  електричному  розряді 
в  газах  є  ефективним  заеобом,  що  застосовується  у  плазмохімії  СОх,  КОх  та 
80х,  використовується  у  багатьох  інших  технологічних  процееах  та 
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наукових  дослідженнях.  У  гірничій  справі  практичну  зацікавленість  і 
важливе  значення  мають  дослідження  причин  утворення  пожежо-, 
вибухонебезпечних  та  агресивних  середовищ  в  підземних  виробках  шахт  і 
рудників  та  розробка  способів  суттєвого  зменшення  впливу  фізико- 
хімічних  реакцій  у  газі. 


1.1.  Хімічні  реакції,  ініційовані  електричними  зарядами, 
і  їх  вплив  на  мікроклімат  усередині  вибухозахищених  оболонок 

Атмосфера  підземних  гірничих  виробок  (температура,  вологість, 
загазованість  і  запилення)  впливає  на  стан,  термін  служби  і  надійність 
вибухобезпечного  електроустаткування.  Але,  додатково  до  зовнішніх 
впливів  гірничої  атмосфери,  в  обмеженому  обсязі  рудникових 
вибухобезпечних  апаратів  формується  особливий  мікроклімат, 
відмінний  по  своїх  параметрах  від  клімату  гірничих  виробок. 
Причинами  цього  є,  наеамперед,  відкрита  електрична  дуга,  що  виникає 
при  комутації  підключеного  до  апарата  навантаження,  а  також  теплові  і 
фізичні  процеси  усередині  вибухонепроникної  оболонки  при 
проходженні  струму  навантаження  по  силовим  і  контрольним  апаратам  і 
струмопроводам,  при  переході  через  окремі  вибухонепроникні 
відділення  і  т.  інш. 

Мікроклімат  у  вибухонепроникних  оболонках  рудникових 
вибухобезпечних  апаратів  керування  і  захисту  характеризується  двома 
групами  факторів:  зовнішніми  і  внутрішніми.  Зовнішні  фактори 
залежать  від  місця  установки  апарата,  тобто  від  температури,  вологості, 
загазованості  і  запилення  гірничої  атмосфери,  а  також  від  характеру  і 
режимів  роботи  підключених  технологічних  машин  і  механізмів.  До 
внутрішніх  факторів,  що  впливають  на  утворення  мікроклімату, 
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відносяться  конструктивні  особливості  вибухобезпечного  апарата  й 
основного  комутаційного  апарата,  тип  і  конструкція  силового 
струмопровода,  спосіб  передачі  електроенергії  з  одного 
вибухонепроникного  відділення  в  інше,  матеріали  і  покриття  контактних 
з'єднань,  тип  і  клас  ізоляційних  матеріалів,  закладених  у  виробі. 

Основні  зовнішні  фактори  і  ступінь  їхнього  впливу  на  мікроклімат 
усередині  вибухобезпечного  апарата  характеризуютьея  такими 
факторами. 

Температура  в  міеці  уетановки  апарата  залежить  від  глибини 
залягання  вугілля,  еистеми  і  етану  провітрювання  гірничих  виробок  і 
екладає  2 5.., З 5  °С.  Вологіеть  атмоефери  зазначених  гірничих  виробок 
складає  відповідно  80...  100%,  запиленість  і  загазованість  гірничих 
виробок  у  місцях  розташування  переважної  більшоеті  рудникових  апаратів 
складає  відповідно  1000  мг/м^  і  1,5  %. 

Вплив  останніх  двох  факторів  при  дотриманні  нормативних 
параметрів  у  гірничих  виробках  еуттєво  не  впливає  на  етан 
вибухобезпечних  апаратів.  Температура  і  вологіеть  навколишнього 
середовища  мають  більш  значний  етупінь  впливу  на  мікроклімат 
усередині  оболонки.  Зв'язок  між  зовнішніми  і  внутрішніми  факторами 
атмоефери  характеризують  наеамперед  її  термодинамічні  параметри, 
оеновними  з  яких  є  об’єм,  тиек,  температура  та  вологіеть. 


1.1.1.  Експериментальні  дослідження 

Температура  повітряної  еуміші  в  оболонці  апарата  змінюється  від 
температури  навколишнього  середовища,  тобто  гірничої  виробки  в  місці 
установки  до  максимальної  сталої  температури  повітря  усередині 
оболонки  при  комутації  підключеного  навантаження  етрумоприймача. 
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Температура  перегріву  повітряної  суміші  усередині  оболонки  апарата 
може  досягати  40...  60  °С.  Нагрівання  мікроатмосфери  усередині 
оболонки  призводить  до  підвищення  тиску  і  витисненню  з 
вибухонепроникної  оболонки  повітряної  суміші,  рис.  1.  Зазначені 
залежності  показують,  що  при  перевищенні  температури  усередині 
оболонки  над  навколишнім  середовищем  на  60  °С  об’єм  витиснутого 
повітря  досягає  10%. 


Рисунок  1  -  Залежність  витиснутого  об'єму  повітря  від  температури 
перегріву  для  оболонок  вільним  об’ємом  (л):  1  -  120;  2  -  240;  З  -  300 


Для  побудови  залежності  відносної  вологості  від  температури 
перегріву  використовуються  експериментальні  дані  з  урахуванням 
табульованих  залежностей  парціального  тиску  від  температури 
середовища.  Криві,  представлені  на  рис.  2,  показують,  що  підвищення 
температури  повітря  усередині  оболонки  апарата  усього  на  10  °С 
призводить  до  зниження  відносної  вологості  майже  в  два  рази.  Ці 
залежності  дозволяють  установити,  що  при  комутації  підключеного 

навантаження  усередині  оболонки  апарата  відносна  вологість  значно 
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нижче  за  вологість  атмосфери  гірничої  виробки.  Однак  при  циклічній 
роботі  апарата  в  короткочасному  чи  повторно -короткочасному  режимах 
роботи  відбувається  періодичне  нагрівання  й  охолодження 
мікроатмосфери  усередині  оболонки,  а  також  зниження  і  підвищення 
відносної  вологості  з  випадінням  конденсату  на  струмоведучих 
елементах  і  ізоляційних  деталях  апаратів. 

Проведені  дослідження  в  лабораторних  умовах  і  в  шахтах 
показали,  що  при  роботі  апарата  в  реальних  режимах  відносна  вологість 
повітряної  суміші  усередині  оболонки  складала  40... 60  %  при  відносній 
вологості  гірничої  виробки  95. ..98  %,  а  перевищення  температури 
мікроатмосфери  апарата  над  температурою  повітряної  суміші  гірничої 
атмосфери  20. ..40  °С. 

Основним  компонентом,  що  визначає  агресивність  мікроклімату  в 
апараті,  є  виділюваний  при  горінні  електричної  дуги  атомарний  азот,  що  у 
з'єднанні  з  киснем  повітря  утворює  окисли  азоту:  0,5-М2+0,5-О2^КО, 
А^«1  еВ. 
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Температура  перегріву,  К 


Рисунок  2  -  Залежність  відносної  вологості  в  оболонці  апарата 
від  температури  перегріву  при  різній  вологості  навколишнього 

середовища 
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Вибухонепроникна  оболонка  апарата  перешкоджає  вільному 
обміну  внутрішньої  мікроатмосфери  з  атмосферою  навколишнього 
середовища.  Окисли  азоту  з'єднуються  з  парами  води,  що  є  присутніми  у 
вологому  повітрі,  і  утворюють  азотну  й  азотисту  кислоти.  Ці  активні 
окислювачі,  осідаючи  на  контактні  з'єднання,  ізоляційні  матеріали, 
металеві  вузли  і  деталі,  приводять  до  окислювання  контактних  і 
рухливих  з'єднань,  зниження  опору  ізоляції.  Збільшення  перехідного 
опору  контактних  з'єднань,  зниження  опору  ізоляції,  затирання  рухомих 
частин  призводять  до  виникнення  відмовлень,  і  виходу  з  ладу  апарата. 

Спочатку  як  комутаційні  апарати,  використовувалися  серійні 
контактори  загальнопромислового  застосування  типу  КТ  6000,  КТ  7000  з 
посиленою  ізоляцією  і  конструктивними  змінами.  Однак  надійність  і 
довговічність  таких  апаратів  виявилася  дуже  низькою  -  термін  їх  служби 
не  перевищував  1...3  місяці.  Для  встановлення  причин  виходу  з  ладу 
апаратів  були  проведені  лабораторні  і  промислові  дослідження,  у  ході 
яких  з'ясувалося,  що  в  результаті  частих  комутацій  підключених 
струмоприймачів  за  рахунок  виділюваних  окислів  азоту  відбувається 
інтенсивний  знос  контактних  з'єднань,  ізоляції  і  вихід  з  ладу  апарата  в 
цілому.  Була  розроблена  нова  серія  спеціальних,  так  званих  вугільних 
контакторів,  що  мають  позначення  КТУ,  у  яких  вперше  у  світовій 
практиці  були  створені  герметизовані  дугогасні  камери. 

Порівняльні  результати  випробувань  для  визначення  залежності 
концентрації  окислів  в  апаратах  з  різними  контакторами  з 
герметизованими  -  КТУ  (крива  2),  негерметизованими  -  КТ  6000  (крива 
1)  і  КТ  7000  (криві  3-4)  дугогасними  камерами  приведені  на  рис.  3. 

При  однакових  параметрах  навантаження  спеціальні  вугільні 
контактори  КТУ  мали  майже  в  три  рази  нижче  концентрацію  окислів 
азоту,  чим  загальпромислові  КТ  при  однаковому  наробітку.  Різна 
концентрація  окислів  азоту  для  контактора  КТ  7000  (крива  З  стосовно 
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кривої  4)  викликана  різним  числом  циклів  "вмикання-відключення"  (В- 
В)  у  годину  (відповідно  600  і  300  В-В  у  годину).  Приведені  на  рис  З 
залежноеті  для  різних  рудникових  пуекачів  евідчать  про  значний  вплив 
на  концентрацію  продуктів  дугогасіння  потужності  підключеного 
навантаження,  чаетоти  комутацій  і  типу  сиетеми  дугогасіння.  При 
використанні  вакуумних  контакторів  сліди  окислів  азоту  в  оболонці 
відсутні 


Кількість  циклів 
"вмикання-відключення"  N-10  ^ 


Рисунок  з  -  Концентрація  окислів  азоту  у  функції  кількості 
напрацьованих  циклів  В-В  для  рудникових  пускачів  ПВИ-320  (крива  1), 
ПВИ-250  (крива  2)  і  для  пускача  ПВИ-125  при  частоті  600  і  300  циклів 

В-В  у  годину  (відповідно  криві  З  і  4) 


Дослідження  показали,  що  виділювана  при  комутації  кількіеть 
окиелів  азоту,  визначена  шляхом  хімічного  аналізу,  значно  нижче,  ніж 
можна  було  б  очікувати,  приймаючи  до  уваги  навантаження,  що 
комутується.  Вихід  окислів  азоту  пропорційний  і  енергії,  що  розсіюєтьея 
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в  дузі  кожної  з  фаз.  Причому  при  включенні  навантаження  виділяється 
менше  1%  окислів  у  порівнянні  з  кількістю  окислів  азоту,  виділюваних 
при  відключенні,  що  цілком  зрозуміло.  Тому  при  дослідженні  виділення 
окислів  азоту  достатньо  враховувати  енергію,  що  розсіюється  в 
електричній  дузі  тільки  при  відключенні  навантаження. 

Виділювана  кількість  окислів  азоту  пропорційна  роботі  (енергії) 
дуги  при  постійному  коефіцієнті  потужності  ланцюга,  а  при  зміні 
коефіцієнта  потужності  змінюється  коефіцієнт  пропорційності  між 
енергією  (роботою)  електричної  дуги  і  кількістю  виділюваних  окислів 
азоту.  На  рис.  4  приведені  залежності  роботи  електричної  дуги  і 
кількості  виділюваних  окислів  від  коефіцієнта  потужності.  Залежність 
коефіцієнта  пропорційності  між  цими  величинами  від  коефіцієнта 
потужності  електричного  навантаження  наведена  на  рис.  5. 


Коефіцієнт  потужності,  С08ф 


Рисунок  4  -  Робота  (крива  1)  і  кількість  окислів  азоту, 
що  виділяються,  (крива  2)  у  функції  коефіцієнта  потужності 


Отримані  залежності  показують,  що  для  рудникових  пускачів  при 

коефіцієнті  потужності  0,3 5... 0,60  (що  відповідає  початковому  періоду 
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пуску  АД),  коефіцієнти  пропорційності  між  кількістю  виділених  окислів 
азоту  і  роботою  електричної  дуги  складає  (7, 5. ..8)  10  ^  г/Дж.  Ці  дані 
отримані  при  роботі  рудникового  пускача  ПВИ-250  у  повторно- 
короткочасному  режимі  категорії  застосування  АС-4.  Значення 
параметрів  роботи  дуги  і  виходу  окислів  азоту  отримані  як  середні 
арифметичні  значення  100  циклів  "відключення". 

Запропонований  коефіцієнт  пропорційності  {рис.  5)  може  бути 
використаний  для  визначення  зміни  концентрації  окиелів  азоту  від 
кількості  напрацьованих  циклів  на  прямолінійній  ділянці  цієї 
характеристики. 


Коефіцієнт  потужності,  С08ф 


Рисунок  5  -  Коефіцієнт  пропорційності  між  роботою  (енергією) 
і  кількістю  виділюваних  окислів  азоту  К^/а  =  8/А  (г/Дж)  у  функції 
коефіцієнта  потужності  навантаження  (електричного  ланцюга) 


Таким  чином,  основним  джерелом  появи  агресивних  компонентів  у 
повітряному  середовищі  усередині  оболонки  є  процее  комутації  (горіння 
дуги  в  повітряному  середовищі),  а  кількість  оксидів  азоту,  що 
утворюється  при  відключенні  навантаження  потужності,  пропорційна 
навантаженню,  залежить  від  її  коефіцієнту  потужності  і  типу  системи 
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дугогасшня,  яка  визначає  інтенсивність  дугогасшня,  та  час  горіння 
електричної  дуги. 

Експериментальні  дослідження  рудникових  вибухобезпечних 
пускачів  показали,  що  мікроюатмосфера  усередині  вибухонепроникної 
оболонки  впливає  на  стан  контактних  з'єднань  контрольних  ланцюгів  і 
ланцюгів  керування.  Під  дією  агресивного  середовища,  що  виникає  при 
комутації  навантаження  контактором  рудникового  пускача,  у 
відділенні  контактора  спостерігається  збільшення  електричного  опору 
контактних  з'єднань,  а  в  деяких  випадках  і  повна  втрата  контакту. 
Оскільки  схема  керування  пускача  забезпечує  контроль  опору  ланцюга 
заземлення  і  ланцюгів  керування,  то  збільшення  опору  контактних 
з'єднань  викликає  помилкові  спрацьовування  захисту  чи  його 
неспрацьовування  в  аварійній  ситуації.  Дослідження  впливу 
мікроатмосфери  на  перехідний  опір  контактних  з'єднань,  виконаних  з 
різних  матеріалів,  дозволяє  прогнозувати  довговічність  і 
безвідмовність  роботи  схеми  керування  і  пускача  в  цілому. 

Для  досліджень  обрані  електричні  штепсельні  контактні  з'єднання 
блоків  керування  і  захистів,  що  були  розташовані  у  вибухонепроникній 
оболонці  рудникового  пускача  ПВИ-125.  Агресивне  середовище 
усередині  оболонки  створювалося  комутацією  струму  350  А  при 
напрузі  660  В  и  коефіцієнті  потужності  навантаження  0,35  з  частотою 
180  циклів  вмикання-відключення  в  годину,  що  відповідало  повторно - 
короткочасному  режиму  категорії  застосування  АС -4.  Критерієм 
відмовлення  групи  контактних  з'єднань  із  трьох  послідовно  з'єднаних 
контактів  було  збільшення  перехідного  опору  понад  10  Ом. 

Візуальний  огляд  і  оцінка  загального  стану  контактних  поверхонь 
після  випробувань  показали,  що  поверхня  контактів  неоднорідна  як  по 
кольору,  так  і  по  характеру  і  виду  скупчень  продуктів  корозії,  рис.  6.  Усі 
покриття  втратили  металевий  блиск,  потьмяніли.  При  п'ятикратному  і 
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десятикратному  збільшенні  виявлено,  що  на  всіх  зразках,  за  винятком 
контактів  зі  срібним  покриттям,  у  місці  контакту  щічок  гнізда  з  ножем 
покриття  частково  зняте  до  латуні. 


Рисунок  6  -  Характерний  вид  контактів  після  випробувань: 
а)  при  п  ’ятикратному  збільшенні;  б)  при  тридцятикратному  збільшенні 


На  всіх  контактах,  у  покритті  яких  міститься  олово,  утворився 
кристалічний  наліт  білого  кольору.  Тому  що  8п02  і  8пО  -  продукти 
корозії  олова,  причому  8пО  на  повітрі  окислюється  до  8п02 ;  на 
контактах  видно  кристали  двоокису  олова  (білого  кольору)  8п02  [19]. 

На  контактах  без  покриття,  або  з  тонким  гальванічним  покриттям 
утворюється  і  сольовий  наліт  бірюзового  або  темно -блакитного  кольору 
у  вигляді  великих  кристалів  взаємодії  пар,  що  є  продуктом  реакції 
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азотної  кислоти  з  міддю,  що  міститься  в  латуні:  Си(  N0^  )2  (білі,  ледве 
зеленуваті  криетали)  і  Си(  N0^  )2  ■  ЗН2О  (темно-блакитні  криетали)  На 
нікельованих  контактах  можуть  утворюватиея  продукти  корозії  N10  і 
М(  N0^  )2  (криеталічні  речовини  зеленуватого  кольору). 

Аналіз  результатів  доеліджень  дозволяє  вказати,  що  протягом  7,5 
місяців  впливу  корозії,  яка  містить  окисли  азоту,  на  поверхні  контактів 
пітепеельних  роз’ємів  з  латуні  без  покриття  і  з  одно-  і  двошаровими 
покриттями  з  нікелю,  або  олова  зі  еплаву  Ві-8п  утворилиея  не 
електропровідні  плівки  продуктів  корозії  у  виді  окислів  і  азотнокислих 
солей  металів  покрить  і  підкладки  (латуні).  Посріблені  ж  контакти 
кородували  порівняно  мало. 

Проведений  хімічний  аналіз  продуктів  корозії  в  оболонці 
рудникового  пускача  після  відмовлення  веіх  його  контактних  з'єднань 
показав,  що  основні  складові  продуктів  корозії  є  окисли  заліза  (66,2%), 
волога  12,3%,  продукти  органічного  походження  (корозія 
лакофарбового  покриття)  21,5%,  маються  сліди  міді  й  олова. 
Інтенсивний  процес  корозії  в  початковий  період  експерименту 
пояенюєтьея  руйнуванням  металу  і  захиеного  покриття  під  дією  окиелів 
азоту  і  їхніх  з'єднань  з  вологою  повітря. 

Утворені  окисні  плівки  мають  високу  стійкість  до  пробою  при 
підвищеній  напрузі.  Напруга  пробою  цих  окиених  плівок  (зниження  опору 
нижче  1  Ом)  екладає  20.. .35  В  для  контактів,  покритих  нікелем  (Л)  і  без 
покриття,  і  40. ..80  В  -  для  інших  покрить.  Окисна  плівка  на  поверхні 
контактів  з  покриттям  нікель  (Л),  олово-нікель  (Г)  не  дозволила  знизити 
опір  контактного  з'єднання  нижче  10  Ом. 

Викориетовуючи  результати  доеліджень  для  опиеу  залежноеті 
перехідного  опору  від  кількоеті  напрацьованих  циклів  була  прийнята 
логарифмічна  функція.  Для  апроксимації  отриманих  даних 
використовується  така  залежність; 
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у  =  а  +  Ьд-  ехр(  Ь^Х  ) , 


де  У  -  перехідний  опір  контактів,  м-Ом;  X  -  кількість  відпрацьованих 
циклів  В-В;  а,  Ьо,  Ьі  -  коефіцієнти  регресії. 

Рудникові  вибухобезпечні  пускачі  в  силу  специфічних  умов 
експлуатації  і  конструкції  істотно  відрізняються  від  електричних 
апаратів  загальнопромислового  виконання.  Вибухонепроникна  оболонка 
рудникового  пускача  утруднює  теплообмін  з  навколишнім 

середовищем,  а  також  викликає  додаткові  втрати  при  протіканні  струму 
по  силових  струмоведучих  шпильках,  що  проходять  через  сталеві 
перегородки  між  вибухонепроникними  відділеннями  однієї  оболонки  й 
інші  конструктивні  елементи,  виконані  з  феромагнітних  матеріалів. 
Дослідження  нагрівання  рудникових  пускачів  [20]  показали,  що  перегрів 
сталевої  перегородки  між  відділеннями  контактора  і  роз'єднувача  в  місці 
розташування  прохідних  шпильок  досягає  95°С  при  струмі  тривалого 
режиму  400  А.  Це  на  30%  перевищує  температуру  повітря  в  оболонці 
рудникового  пускача  в  зоні  розташування  прохідних  шпильок. 

Використовуючи  експериментальні  дані  по  нагріванню 

рудникових  пускачів  приймаємо  наступні  допущення:  температуру 
навколишнього  повітря  усередині  вибухобезпечної  оболонки  пускача 
приймаємо  рівною  353  К  (80°С),  температуру  навколишнього 
середовища  303  К  (30°С),  ступінь  чорності  випромінювання  для  сталі 
і  чавуна  є  =  0,80.  За  визначальний  розмір  для  шайб  і  інших  подібних 
елементів  приймаємо  їхню  товщину,  для  перегородок  між 
вибухонепроникними  відділеннями  через  їхні  значні  розміри 

визначальний  розмір  не  входить  у  формулу  для  визначення 

сумарного  коефіцієнта  тепловідводу,  що  визначається  по  наступній 
залежності: 
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де  Я  -  коефіцієнт  теплопровідності,  Вт-м  ^трад 

середня  температура,  К. 

Перевищення  температури  окремих  елементів  рудникового 
пускача  можна  представити  в  загальному  виді  таким  виразом: 
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де  магнітна  проникність  для  конструкційної  сталі;  /=50  - 

частота  мережі,  Гц. 

Для  конкретного  елемента  конструкції  усі  величини  крім 
струму  і  магнітної  проникності,  можна  замінити  деякою  постійною  і 
представити  залежність  температури  перегріву  у  виді  більш  простої 
залежності: 


АТ  = 
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де 


А  = 


71  -  у  -  к^  • 


Аналіз  показує,  що  розбіжність  розрахункових  та 
експериментальних  даних  становить  не  більше  12%,  табл.  1. 
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1.1.2.  Утворення  оксидів  азоту  у  нерівноважній  плазмі 


Швидкість  коливального  збудження  молекул  азоту  є  достатньо 
високою  і  тому  основна  доля  розрядного  енергетичного  внеску 
локалізується  власне  на  коливальних  ступенях  евободи,  що  забезпечує 
виеоку  енергетичну  ефективніеть  процееу.  Для  таких  молекул  як  СО2, 
N2,  Н2О  при  температурі  електронів  до  З  еВ  оеновна  доля 
зосереджується  на  збудженні  коливальних  ступенів  свободи  основних 
електронних  етанів.  Ветановлено,  що  коливальна  енергія  найбільш 
ефективна  у  подоланні  активаційних  бар'єрів  ендоергічних  реакцій. 

Окислення  повітряного  (молекулярного)  азоту  у  плазмі  є  одним  із 
найвідоміших  плазмохімічних  процесів  тому,  що  молекули  азоту 
інтенеивно  збуджуютьея  коливально  під  час  електронного  удару  і 
характеризуютьея  дуже  низькою  швидкіетю  коливальної  релаксації  [15; 
16].  Утворення  в  плазмі  оксидів  азоту  є  типова  ендоергічна  ланцюжкова 
реакція,  що  етимулюєтьея  коливальним  збудженням  основного 
електронного  етану  молекул  в  плазмі. 

Таблиця  1 

Параметри  апроксимації  та  перегріву  елементів  конструкції 
рудникового  пускача  ПВІ-320 


Наймену¬ 

вання 

елемента 

Коефіцієнти 

апроксимації 

Перегрів 

розрахунковий  /  експериментальний, 

°С,  при  струмах  навантаження  у  А 

Аі, 

А^-Гн  / 
(м-град^) 

ао, 

Гн/м 

аь 

Гн/  (м-А) 

100 

200 

300 

400 

Похибка, 
тах  % 

Сталева 

перегородка 

1126 

ІМО'^ 

2-10'’ 

8,7/9 

31,5/35 

59,9/60 

86,2/90 

11Д 

Центрувальна 

шайба 

299 

8,67- 10'^ 

1,67- 10'^ 

15,3/16 

53,4/48 

99,5/90 

130,5/135 

11,2 
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Утворення  N0  у  плазмі  визначається  набором  елементарних 
процесів,  серед  яких  лімітованими  є  такі,  де  проходе  руйнування 
міцного  зв'язку  (~10  еВ)  молекули  N2.  Механізм  руйнування  зв'язку 
досліджувався  у  даній  роботі  з  використанням  квантово-механічної 
задачі  про  плив  електрона  на  два  кулонівські  центри.  Коефіцієнт 
корисної  дії,  що  був  одержаний  в  експериментах  по  синтезу  N0  складає 
4... 20%,  а  енерговитрати  -  (4... 25)  еВ/моль.  Енергетично  найбільш 
ефективні  реакції  дисоціації  молекул  азоту  (як  і  СО2)  зв'язані  з 
коливальним  збудженням  азоту,  де  концентрується  основна  частка 
енергетичного  вкладу  в  розряд  [15]. 

Як  відомо  [16],  елементарна  реакція  синтезу  N0  за  участю 
коливально-збуджених  молекул  азоту  лімітується  елементарним 
процесом 

0(^Р)  -к  К2*(^ 2/’^)  ^  N0(^71)  -к  N(^8),  Еа^ЗеВ, 

що  визначає  особову  зацікавленість  до  його  аналізу.  Приведена  реакція 
може  здійснюватися  як  прямим  адіабатичним  так  і  неадіабатичним 
шляхом  [15].  Імовірність  адіабатичного  шляху  одержують  із  принципу 
детальної  рівноваги.  Неадіабатичний  шлях  синтезу  здійснюється  через 
утворення  проміжного  коливально-збудженого  комплексу  М02*(^22^): 

0(^Р)  -к  М2*(^5:2^у)  ^  ^  N0(^71)  -к  N(^8)  (1) 

Рекомбінація  (перша  стадія  процесу)  характеризується  лише 
імовірністю. 

Крайнє  значення  енергетичної  ефективності  утворення  N0  у 
нерівноважних  розрядах  не  перевищує  35%.  У  звичайних  умовах  роботи 
апаратів,  що  відрізняються  від  експериментальних  умов,  ефективність 
синтезу  наближується  до  10%. 
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друга  стадія  -  руйнування  проміжного  комплексу  N02*(^I1/)  - 
описується  закономірностями  статистичної  теорії.  У  рамках  класичного 
наближення  Маркуса-Райеа  з  урахуванням  поправок  на  ангармонізм, 
імовірніеть  руйнування  зв'язку  буде  залежати  від  коливальної, 
поступальної  і  обертальної  енергії.  Розрахунки  [15]  показують,  що  на 
швидкість  реакції  впливає  переважно  коливальна  енергія. 

Таким  чином,  в  умовах  нерівнованості  плазми  при  температурах  до 
1000  К  може  переважно  здійснюватися  механізм  (1)  елементарного 
процесу  синтезу,  який  принципово  відрізняється  від  такого  в  рівноважних 
умовах  [15;  16].  Але,  якщо  у  термодинамічно  рівноважній  сиетемі  будуть 
знаходитися  електричні  заряди,  то  і  в  останньому  випадку  без  зовнішньої 
енергетичної  дії  елементарний  процес  синтезу  буде  йти  шляхом  (1). 
Тобто,  іншими  словами,  молекули  азоту  (і  навіть  оксиду  вуглецю) 
можуть  дисоціювати,  якщо  на  хімічний  зв'язок  буде  здійснюватися  вплив 
зовнішнього  поля  точкового  заряду  (іону)  -  кулонівського  центру. 

Іншим  можливим  джерелом  енергетичного  впливу  на  ініціювання 
хімічних  реакцій  може  бути  післясвічення  (зона  релакеації  параметрів 
плазми),  що  виникає  безпосередньо  за  електричним  розрядом  у  повітрі. 
Спектральними  дослідженнями  показано  [21],  що  впливи  розряду 
супроводжуються  не  тільки  ерозією  електродів,  нагріванням  та  іонізацією 
молекулярного  газу,  але  й  збагаченням  екладу  газу  хімічними 
компонентами  з  високою  ентальпією  утворення  і  збудженим  станом 
атомів  і  молекул.  Параметри  післясвічення  розряду:  час  ~  0,1  мс,  газова 
температура  ~  2000. .  .3000  К. 

Результати  досліджень  [22;  23]  показують,  що  характерний  час 
спадання  іонізації  у  початковій  стадії  післясвічення  корелює  з  часом 
розпаду  метастабільних  станів.  У  роботі  [22]  встановлені  параметри 
залежності  швидкості  іонізації  у  післясвіченні  розряду.  Автори  роботи 
роблять  висновок  про  те,  що  стадією  процесу  іонізації  в  самостійному 
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розряді  в  азоті  в  умовах  сильної  коливальної  нерівнованності  є  зіткнення 
молекули  у  високозбудженому  етані  з  електронно-збудженою  молекулою 
в  нижчележачому  триплетному  стані . 

Додавання  N02  в  зону  післясвічення  [24]  супроводжується 
люмінесценцією  у  зелено-блакитній  частині  спектру,  яка  пов’язана  з 
океидом  міді  (матеріал  електроду),  і  пояснюється  як  результат 
хемілюмінеецентної  реакції  між  атомами  міді  і  молекулами  двоокису 
азоту.  Додавання  в  зону  післясвічення  кисню  приводило  до  суттєвого 
збільшення  емісійного  випромінювання  атомів  міді,  пояснене  [25] 
приекоренням  передачі  енергії  від  збуджених  молекул  азоту  атомам  міді. 

Триваліше  післяевічення  виникає  в  обмеженому  міжелектродним 
проміжком  об'ємі  молекулярного  газу,  у  тому  числі  повітря,  азоту,  кисню 
або  закису  азоту,  в  результаті  короткочасної  дії  потужного 
високочастотного  розряду,  для  якого  характерне  нерівноважне 
енерговкладання  [26]. 

Аналізуючи  вищезазначені  роботи,  логічно  припустити,  що 
негативний  вплив  дугогасильних  процесів  у  автономній  атмосфері 
вибухозахищених  оболонок  може  поеилюватися,  або  зменшуватиея  в 
залежності  від  ступеню  відхилення  тих  чи  інших  компонентів  газу  від 
екладу  звичайної  повітряної  атмосфери  або  наявності  домішок,  що 
можуть  бути  потенційними  каталізаторами  хімічних  реакцій.  Це,  так  би 
мовити,  є  додаткові  фактори,  вплив  яких  може  бути  потенційно 
негативним,  з  точки  зору  утворення  нових  токеичних  і  агресивних 
з'єднань. 

Таким  чином,  вирішення  проблеми  надійноеті  роботи  шахтних  та 
рудникових  апаратів  захиету  та  управління  технологічними 
комплексами  слід  вирішувати  з  урахуванням  результатів  дослідження 
фізико-хімічних  ефектів. 
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2.  ФІЗИКО-МАТЕМАТИЧНА  МОДЕЛЬ  РОЗРАХУНКУ 
ЕЛЕКТРОННИХ  ТЕРМІВ  МАЛИХ  МОЛЕКУЛ 


В  основу  рішення  задачі  взаємодії  компонентів  молекулярних 
систем  з  точковими  електричними  зарядами  покладена  квантово- 
механічна  модель  визначення  стану  електрона  в  полі  двох  кулонівських 
центрів  [27-29].  Викориетовується  адіабатичне  наближення  Борна- 
Оппенгеймера  (рух  електронів  розглядаєтьея  окремо  від  руху  ядер, 
відстані  між  якими  вважаються  фіксованими). 

Для  вирішення  двоцентрової  задачі  з  рівняння  Шредінгера  виділено 
модельне  рівняння,  рішення  якого  одержано  в  аналітичному  виді.  Облік 
граничних  умов  дозволив  знайти  квантові  числа,  що  визначають  стан 
валентного  електрона  в  полі  двох  кулонівських  центрів.  Одержане  рішення 
є  наближеним,  оскільки  воно  представляє  рішення  укороченого  рівняння 
Шредінгера  і  тому  може  вважатися  нульовим  наближенням  до  точного 
рішення. 

Наступним  логічним  кроком  був  облік  невживаних  у  модельному 
рівнянні  членів  точного  рівняння  шредінгера  як  збурення  до  нульового 
наближення.  Вибраний  підхід  реалізований  за  допомогою  знайденої  нами 
функції  Гріна  [28],  що  є  рішенням  модельного  рівняння  з  дельтаподібними 
джерелами.  Знання  функції  Гріна  дозволило  визначити  в  аналітичному 
вигляді  хвильову  функцію  оеновного  стану  в  першому  наближенні.  Пошук 
наступних  наближень  в  аналітичному  вигляді  недоцільний  з  двох  причин: 
рішення  дуже  громіздкі  і  не  дають  іетотної  поправки  до  першого 
наближення. 

Кулонівський  потенціал,  що  впливає  на  хімічний  зв'язок  молекули, 
розглядається  як  збурення  до  модельного  рівняння  шредінгера.  Тому 
була  застосована  функція  Гріна  для  знаходження  аналітичного  рішення, 
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що  було  одержано  в  еліпсоїдній  системі  координат,  в  якій  зовнішній 
потенціал  представлявся  у  вигляді  розкладання  Неймана. 

Для  визначення  параметрів  морсівського  потенціалу  і  вирішення 
коливальної  квантово-механічної  задачі  вперше  розроблено  підхід, 
заснований  на  апроксимації  електронного  терму  потенціалом  Морса. 

Використовування  рішень  квантової  механіки  для  знаходження 
макропараметрів  досліджуваних  систем  засновано  на  вживанні  методів 
статистичної  фізики,  а  саме,  -  розподілу  Гіббса.  В  цьому  випадку 
енергія  молекулярної  системи  представляється  у  вигляді  суми 
електронного,  коливального,  обертального  і  поступального  рухів.  Із 
статистичної  суми  знаходяться  розподіли  електронів  по  енергетичних 
станах  залежно  від  температури  системи.  Іншими  словами, 
статистичний  підхід  встановлює  взаємозв'язок  квантово-механічних 
рішень  з  макроскопічними  параметрами  системи. 

В  загальному  випадку  в  стаціонарному  рівнянні  Шредінгера 
повний  нерелятивістський  гамільтоніан  молекули  можна  представити  у 
вигляді  суми  трьох  операторів; 

Н  =  -к  Т\г)  +  и(г,К),  ( І ) 

де  і  -  кінетична  енергія  ядер  і  електронів,  ІІ(г,К)  - 

потенційна  енергія  молекули,  що  включає  енергію  тяжіння  електронів 
до  ядер,  енергію  міжелектронного  та  між'ядерного  відштовхування,  К  і  г 
-  ядерні  і  електронні  координати. 

Строго  кажучи,  рухи  електронів  і  ядер  у  молекулі  зв'язані  між 
собою,  проте,  враховуючи  той  факт,  що  маси  ядер  в  1 0^-10^  разів 
більше  маси  електронів  і,  отже,  швидкості  руху  ядер  у  багато  разів 
менше  швидкостей  електронів,  подальші  розрахунки  проводитимуться  в 
адіабатичному  наближенні  шляхом  розділення  цих  внутрішньо 
молекулярних  рухів,  тобто,  припускаючи  їх  незалежний  характер.  Тоді 


ЗО 


підсистеми  електронів  1  ядер  характеризуватимуться  кожна  своїм 
рівнянням  Шредінгера  і  своєю  хвильовою  функцією. 

Необхідно  відзначити,  що  точне  рішення  рівняння  шредінгера 
вдається  одержати  тільки  для  сиетеми,  що  складається  з  двох  частинок  - 
протона  та  електрона  (атом  водню).  Для  розрахунку  молекулярних  термів 
в  даний  час  розроблено  багато  підходів,  найпопулярнішим  з  них  є  метод 
лінійної  комбінації  атомних  орбіталей  (ЛКАО)  [30-33].  Розрахунки 
молекул  Н2,  О2,  N2  і  т.  д.,  указують  на  задовільну  подібність  з 

результатами  експериментів.  Проте  використовування  даного  методу,  як  і 
інших,  вимагає  проведення  великих  чисельних  розрахунків,  аналіз  яких 
утруднений  у  зв'язку  з  відеутністю  аналітичних  рішень.  Тому  в  даній 
роботі  розробляється  модель  розрахунку  молекулярних  термів,  заенована 
на  використовуванні  рішень  задачі  про  стан  електрона  в  полі  двох 
кулонівських  центрів.  За  базис  приймаються  двоцентрові  хвильові 
функції  в  еліпсоїдних  координатах,  всі  функціонали  обчислюються 
аналітично,  а  побудовані  замкнуті  вирази  для  функцій  Гріна  дозволяють 
вирішувати  квантово-механічну  задачу  зі  збуреннями. 

Відомий  метод  ЛКАО  не  дає  можливості  проводити  подібні 
розрахунки  тому,  що  в  ньому  хвильова  функція  визначається  виразом 
\\іі  =  ЕСі^•ф^(і)  (Су  -  коефіцієнти,  що  розраховуються  в  кожній  точці 
простору).  Тут  же  виникають  труднощі  з  урахуванням  енергії  електрон- 
електронної  взаємодії.  Таким  чином,  при  рішенні  задачі  про  стійкість 
реальної  структури  речовини  слід  розвивати  розрахункову  модель, 
запропоновану  в  [27-29;  36]. 

Необхідно  відпрацювати  схему  для  розрахунку  електронних  термів 
будь-яких  двоцентрових  молекул.  Якщо  представити  двоатомну 
молекулу  як  два  позитивно  заряджені  кулонівські  центри  в  оболонці  з 
електронів,  які  кожний  з  атомів  віддав  на  утворення  хімічного  зв'язку,  то 
нашу  сиетему  зарядів  можна  зобразити  у  вигляді  деякого  "кокона"  з 
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еліптичною  симетрією,  при  цьому  хвильова  функція  валентного  електрона 
описує  його  стан  в  полі,  створюваному  зв'язаними  електронами  і  зарядами 
атомних  ядер,  що  створюють  молекулярний  остов.  Вводиться  ефективний 
потенціал,  під  яким  мається  на  увазі  ядерний  потенціал,  екранований 
електронним  остовом  зв'язаних  електронів.  При  розрахунках  електрон- 
електронної  взаємодії  використовується  слеттерівська  константа 
екранування,  яка  враховує  той  факт,  що  заряди  зв'язаних  електронів  не  є 
точковими,  а  "розмазані"  в  просторі  поблизу  ядра.  У  разі  великих  остовів 
зв'язаних  електронів  вони  поляризуватимуться  за  рахунок  взаємодії  один  з 
одним  на  малих  між'ядерних  відстанях.  Цю  поляризацію  також  необхідно 
враховувати  при  проведенні  розрахунків. 


2.1.  Аналітичне  рішення  модельної  задачі 


Розглядається  система  трьох  частинок  із  зарядами  2 1,  22,  2з  і 
масами  Мі,  Мг,  Мз,  взаємодіючих  за  законом  Кулона.  Тут  заряди  2і, 
22  >  0,  а  2з  <  0.  Гамільтоніан  для  електрона  в  полі  двох  кулонівських 
центрів,  при  введенні  координат  Якобі,  запишеться  таким  чином: 


ТЛ-  ^  А  ^  А 

Н  — - Аш - 

2Мо  2М 


2і  •  22  2і  •  2з  2з  •  22 

- А^+— ^ - - - - - ^ - ^ (2) 

2т  К  Га 


де  Мо  =  Мі  +  М2  +  Мз;  А  =  +  ±  =  + ^ _ ; 

М  Мі  М2  ’  т  Мз  М1+М2’ 

г  =  ІК.  —  К.  1 1 

^  “  відстані  від  іонів  а  і  Ь  до  даної  точки  простору 

=  |Кз  —  К2 1] 

валентного  електрона,  К  -  між'ядерна  відстань. 

Вважатимемо,  що  Мз « (Мц  М2),  тобто  частинка  з  масою  Мз  і 
зарядом  2з  може  бути  електроном  або  мезоном.  Надалі  розглядається 
задача  за  участю  електрона.  Перший  член  гамільтоніану  (2)  описує  рух 
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центру  мас  системи  трьох  частинок,  другий  -  відносини  рух  частинок 
Мі  та  М2. 

В  атомній  системі  одиниць  /2  =  23  =  Мз  =  1  та  в  еліпсоїдній  системі 
координат  рівняння  Шредінгера  для  задачі  про  рух  електрона  в  полі 
двох  нерухомих  ядер,  видалених  на  відстань  К,  має  вигляд: 

а"  І 

(3) 

+  2  •  [Е  +  =  0, 


де  А,  ^  ^  ,  1  <  А,  <  ^;  =  -1<А.<1;  П(Л-,ц) 

К.  К. 


1.(^  +  _52_) 

Я  А-  +  ц  А-  —  ц 


Га  і  Гь  -  відстані  до  електрона  від  ядер  2і  і  22  відповідно;  ІІ(А,,|и)  - 
потенційна  енергія  електрона  в  полі  зарядів  2а  і  2ь;  К.  -  між'ядерна 
відстань  атомів  молекули. 

Рішення  рівняння  (2)  шукається  у  вигляді: 


=  Х(?і)-У(ц)-ехр(іЛф),  0<ф<27г.  (4) 


Оскільки  4^(?і,|іі,ф)  =  4^(А,,|іі,ф  +  27т),  то  умовою  квантування  по 
квантовому  числу  Л  буде:  |Л|  =  0,  1,2,  ...  . 

Якщо  підставити  (2)  в  (3),  то  одержимо  систему  двох  звичайних 
диференціальних  рівнянь  другого  порядку 


д 


(. 


'І: 


Л" 

1-ц2 


ц^є  +  |т2  —  А 


У(ц)  =  0, 


(5) 


а 

Ж 


л" 

-1 


+  Х^&  +  22+ 


+  А 


Х(2)=0, 


(6) 
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де  А  -  постійна  розділення;  є 


=  2М=(2,±2ь).К,  де  2  -  є 

ПИТОМИМ  дипольним  моментом  системи. 

Для  розвитку  аналітичного  формалізму  зручно  провести  наступну 

заміну  змінних:  А.  - 1  =  і,  х(і)  =  [і  •  (і  +  2)]~2  •  і'(і)  де  0  <  І  <  оо.  Тоді  (5) 
можна  перетворити  до  вигляду: 


1-Л^  С 


ді' 


41^  ї 


ґ(0+М0^(0=о- 


(7) 


де  \У(і)  = - 

1  +  2 


Л2-1  А  є  2+ 

- 1 - 1 - 

4  2  2  2 


1-А2  ^  А2-1  А  є  2+ 

+ - С  = - +  — +  -  +  — 

4(1  +  2)"  4  2  2  2 


Заміна  у  =  2  ,/-  є  ■  і  приводить  до  наступного  рівняння 


Г  1  к  1-А^^ 


ду‘ 


+ 


-  —  +  —  + 


4  у  4у^ 


ґ(у)+Л¥(у)ґ(у)  =  0. 


(8) 


Тут  к  =  ^  а  \у(у)^  \¥(і;)  при  вищезгаданій  заміні. 

2л/—  є 


Приведення  рівняння  (6)  до  вигляду  (7)  продиктовано  лише  тим, 
що  укорочене  рівняння: 


Г 


ду‘ 


+ 


1  к  1-А 
-  +  + 


2  Л' 


V 


4  у  4у' 


У 


Ду)  =  0 


(9) 


розв'язується  аналітично,  а  два  його  лінійно  незалежні  рішення  є 
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функціями  Уїттекера  М(у)  д  і  ^V(у)  д  [10].  Оскільки  функція  М(у)  д 

к,—  к,—  к,— 

2  2  2 

регулярна  в  точці  у  =  0,  то  такій  математичній  моделі  відповідає  фізична 

задача,  у  якої  ядра  знаходятьея  в  етійких  положеннях.  Оекільки 

збурення  не  має  особливоетей  в  інтервалі  [0,  оо]  і  Ііт  "НУ  =  0 ,  то 

у^оо 

безперервна  чаетина  спектру  буде  однакова  для  рівнянь  (8)  і  (9). 
Оекільки  нас  цікавлять  рішення,  регулярні  в  нулі  і  на  нескінченноеті,  то 
з  явного  виразу  для  функції 


М(у) 


к, 


Л 


Л+1 

=  У  ^ 


_у 

2 


Л  +  1 


1!(Л  +  1)  2!(Л  +  і)(Л  +  2) 


+ ... 


необхідно  зажадати  виконання  наступних  співвідношень: 


2М  +  1 


Л  =  2т; 


9 


\^к  =  N  Л  =  2т  + 1; 


де  т=0,  1,  2;  т<К.  Таким  чином,  разом  з  цілими  квантовими  числами  є  і 
дробові. 

Розглянемо  випадок,  коли  =  0,  або  2^  =  С  і  Л  =  1 .  Оетанні  дві 

рівності  означають,  що  А  +  і  =  ^  і  ^+^  +  ^  =  2к  -є  або  2^=2кл/-є,  а 

2  2  2  2  2  2 

=4к^(-є).  Тоді  Е  =  ;  де  к=  1,  2,  3.  і,  отже,  задача  без 

2к^ 

урахування  збурень  є  воднеподібною. 

Рівняння  (5)  представимо  у  вигляді: 
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(10) 


^  (і  -  ц"  )^  -  +  А  ¥(ц)  +  \¥(ц)у(ц)  =  О, 

д]л  д]л  І- 

де  \¥(А)  =  АІ^  -  А^ІІ;  \¥(А)  -  безперервна  функція  на  відрізку  [-1,  1]. 
Рішеннями  укороченого  рівняння  (9) 

д  2^  5 

—  (І-ц  ) - 

д\і  д\л  І- уі 

при  А  =  п-(п+1)  є  приєднані  функції  Лежандра  рЛр)  і  ^Ар)•  Але  оекільки 
|и  =  [-1,  1],  то  в  подальших  розрахунках  нами  викориетовуватиметьея 
клас  функцій  {Рі^}. 

Аналіз,  проведений  для  рівнянь  (5)  і  (6),  дозволив  виділити 
модельну  задачу,  якій  відповідає  рівняння 


1- А^  с  а  б  2\д  А^ 

ак^  41^  і  діл 

,  (12) 

+  п(п  +  1)  +  —  +  А  =•  •  Р(і,  ц,  ф)  =  0 

5ф  ] 


Для  повернення  від  математичної  моделі  до  реальної  фізичної 
задачі  перетворимо  рівняння  (12),  шляхом  зворотної  заміни  А.  =  і  +  1, 
тоді  ДА.)  =  [Х^  -  1]"^^^-Х(?і).  Оскільки  це  перетворення  нелінійне, 
одержуємо  наступне  рівняння: 
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4  д^  _  4(к  -  '2)^~\x  _  2ц^Е(к  -  2) 


(13) 


ккк -П1-|л^  Уф^  кУ-|д2 


Х2-ц2 


У  к  ■  У (*  - 1^"  )^ + Ут  У |х(^)^(й)Ф(ф) = о 

К^ІУ-ц^)[Ф  д\х  1-ц2  0ф^]1 


Видно,  що  рівняння  (13)  включає  початкове  рівняння  (3).  Таким 
чином,  рішення  модельної  задачі  приймає  вигляд: 


1  —  ^  Гл  1  Л 

х,/  =  Р  =  [і(і  +  2)]“2 -у  2  .е2.ф  к+  ,Л  =  1,у  РУ(іи)е'^^.  (14) 

V  2  2  ) 


де  Ф  -  вироджена  гіпергеометрична  функція. 

Розглянемо  першу  групу  енергетичних  термів  (при  п  =  0)  для 
випадку,  коли  Л  =  0.  Введемо  а^  =  2Ьі  =  2^-  .  Тоді  з  рівняння  (14), 

одержимо  явний  вид  модельних  хвильових  функцій  УкЛ.п: 


Тут  Ьі  знаходимо  за  допомогою  формули 
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Ьі  =-2к±  ,/4кі +2Ь  (18) 

т  _  - 1  ,ч  2+ 

4  2 

Враховуючи,  що  для  зв'язаних  станів  Ьі  >  0,  у  виразі  (18) 
вибирається  знак 

На  основі  рішень  модельної  задачі  одержуємо  енергії  станів, 
відповідні  квантовим  числам  к  =  5/2,  3/2  та  1/2,  обчислені  по  формулі 


'к,  Л,  п 


,Ук,Л,п  1^0  |Ук,Л,п  / 
V  к,Л,п  V  к,Л,п  / 


І  Я  1  І  ^ 

Е  ^  =2+  16Ь^  +  28ЬІ  +  14Ь^  +  -  + - --4Ь^-  21ЬІ  -  40Ь^  -  16Ь^  +  х 

|,0,0  ^  ^  ^  ^  2  2Ь5  )  8  ^  ^  ^  ^ 


д(-4Ь5)[; 


2(б4Ь5  +12885  +8085+1685  +і)-648^  -17685  -1448^  -З685  -Ь5  +^ 


2  5 


-І-уЬз  -^ь|  -^ьу  5є‘">=Еі(-4Ь5/і6ь|  +32ЬЗ  +20Ь|  +405,  +2 


^.0.0  9 - 4Ьз  _ібЬ^ -4Ьз(і  +  4ЬзЄ  •е'"’ЗЕі(-4Ьз) 

О 


ХІ2.  І  ^  +  Ь|  1-2"Гі  +  Ьз  +4Ьз1-є4^зе.(-  463): 


X  1Ьз-Ь^-4Ь"+2+(і6Ь^+8Ь?+Ьз)|; 
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(21) 


Е 


1 

2 


,0,0 


і  (Ьі  -  2+  )+  Еі  (-  4Ьі  (ь2  -  Ьі2^ 


2Ьі  З  '  IV  1/ 


Якщо  в  одному  і  тому  ж  квантовому  стані  (к,  Л,  п)  знаходяться  два 
електрони,  то  необхідно  враховувати  їх  взаємодію  один  з  одним. 
Хвильова  функція  сиетеми  двох  електронів  може  бути  предетавлена  у 
вигляді  детермінанта  Слетера 


'к  (іеі 


иі(і) 

“і  (2) 


“2(1) 

"2(2) 


(22) 


де  Ні  =  аф(1);  и2  =  Рф(2);  ф(1)  і  ф(2)  -  базиені  функції,  які 
викориетовуютьея  в  якоеті  рішення  модельної  задачі  (15)  -  (17). 
Позначення  (1)  і  (2)  є  еукупніетю  координат  першого  і  другого 
електронів  відповідно.  Для  того,  щоб  хвильова  функція  фііеі  була 
антиеиметричною,  необхідно  поклаети  а  =  р.  Тоді,  енергія  електрон- 
електронної  взаємодії  визначитьея  з  виразу: 


Е(ее)к,л,п 


Леї 


'Т  Леї 


\ 


\ 

/ 


Леї 


(23) 


де  1/Гі2  -  енергія  кулонівеької  взаємодії  електронів,  яку  можна 
представити  у  вигляді  розкладання  Неймана  в  еліпсоїдній  системі 
координат;  Гі,2  -  відетань  між  електронами. 

Рішення  реальної  задачі  має  вигляд: 
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(24) 


Ук,Л,п  ф)  =  Фк,Л,п  ф)  + 

+  І  р{і,  Ґ,  Ц,  ]л\ ф,  11і'ф')Ч^к  Д  „  (і',  ]л\ ф')ат'. 

Збурення,  що  доводить  дану  модель  до  реальної  двоцентрової 
задачі,  записується  тепер  таким  чином 

+  2 

\У(і,  ц)  =  \У(()  +  \У(ц)  =  у  +  ц2- ,  (25) 

Оскільки  збурення  є  безперервна  і  обмежена  функція,  то 
ітераційний  процес  сходитиметься. 

Хвильова  функція  в  першому  наближенні  визначається  по  формулі 


(і,  ц,  ф) = \\/{и  ц,  ф) + 1  о{и  ц,  ц',  ф,  фО  \л\ ф')  \\/{и  ц,  ф)іт. 

Рішення  хвильової  функції  з  квантовими  числами  1/2,  0,  0  має 
такий  вигляд: 


1  і  А 

^  1/2  о  0 ^  “^е^  {а^і  +  А2  1п(і  +  2)  + Аз  +  Х; - 

’  ’  г=і(і  +  2У 


+  2  ц  +  уц' 


В]^1  +  В3  +  X 


В, 


Г=і(і  +  2У  ] 


(26) 


де  Аі  =  -2+  -  4аМе'^“Е_і  (-4а)  +  ^  +  (-4а)(6а^  +  ^) ; 

Аг  =  -22  +  е^“Еі(-4а)  +  1/4; 
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Аз  =  —22+  - - 

а.  а. 


-  2(3  +  4а  -  С)е"^“Еі  (-  4а)  -За  +  2-  -1п2- 


11а  +  12а^  +  -  — ЗаС 
2  4 


(-  4а)  -  22  +  е"^“Еі  (-  4а)  х 


4схп 


X 


ОО  г-1  (— і)^В„  °°  Г  г-1  (_  іУ  +  Р  їв 

-С-ІП2+ Х  Е— ^ Е  Е  Е - ^ ^ - 

г=і^=і(4а)’^  (г  —  іу.  г=-і  ^=ір=і(4а)’^  (г  —  ^). . .  (г  —  ^  —  р  +  і) 


£г^і  (-1)^' 

г=і^=і(4а)’' (г  -  і) 


;  Аг  =  Е 


■■^1  (-іУ'в, 

і-і{2а.У  {г  -  еуЧ 


За 


1 


■2: 


Ві  =  -4ае^“Еі(-4а);  Вг  =  -2ае'^“Еі(-4а); 


Вз  -2  +  2(3  +  4а-С)е'^“Еі(-4а) 


ОО  г-1г-^  (-іУ+Р  г-І  (- іУ  ^ в ^ 

Е  ЕЕ - ^ ^ +  ^ 

г=-і^=ір=і(4а)^^Р^^(г  —  £)...{г  —  £  —  р  +  і)  г=Е=і(4а)^(г  -  £) 


ОО  г 


-1  (-іУ"^Вг 

-С-ІП2  +  X  ^ 

г=Е=і(4а)^  (г  —  £)і 


_г-і  (-і)''Вг 

^5і(2а)''(г-^)2^  ' 

Розрахунки  показують,  що  знаходження  хвильових  функцій  більш 
високого  порядку  недоцільне,  оскільки  виходять  дуже  громіздкі  вирази, 
а  поправки  до  енергій  електронних  термін  є  незначними. 


2.2.  Екранування  ядерних  потенціалів  зв'язаними  електронами 

Квантово-механічна  модель  представляє  двоатомну  молекулу  як 
два  позитивні  кулонівські  центри  із  зарядами  2а  і  2ь,  що  зв'язані 
валентними  електронами.  Таке  наближення  не  враховує  наявність 
електронів,  що  екранують  потенціали  ядер.  Для  спрощення  розрахунків 
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молекулярних  систем  можна  вважати,  що  внутрішні  електрони,  які  не 
беруть  участі  в  утворенні  хімічного  зв'язку,  знаходяться  в  станах,  що 
характерні  для  ізольованих  атомів  і  їх  роль  зводиться  лише  до 
екранування  потенціалу  ядра.  Іншими  словами,  електричний  потенціал, 
створюваний  окремим  іоном  молекули,  можна  представити  у  вигляді 
суми  потенціалів: 

2 

ф  =  — +  Фе,  (27) 

г 

де  -  фе  сумарний  потенціал,  що  створюється  2-2у  зв'язаними 
електронами;  2  -  заряд  ядра  (кількість  протонів  ядра);  2у  -  число,  рівне 
кількості  електронів,  що  залишили  атом  для  утворення  хімічного  зв'язку 
молекули. 

природно  очікувати,  що  на  великих  відстанях  від  ядра  фе^О,  а  на 
малих  відстанях  фе  приймає  кінцеве  значення.  Необхідність  введення 
потенціалу  фе  обумовлена  тим,  що  позитивний  заряд  протонів 
локалізований  в  малій  області  ядра,  а  заряд  внутрішніх  електронів  атома 
розподілений  поблизу  ядра  в  області,  розміри  якої  значно  перевищують 
розмір  ядра. Потенціал  фе  визначається  з  рівняння  Пуасона 


Афе  =  -  4пр,  (28) 

І  1 2  І  ... 

де  р  =  -  X  Фі  -  густина  заряду  зв  язаних  електронів;  фі  -  їх  ХВИЛЬОВІ 
і=1 

функції. 

Якщо  функцію  фе  представити  у  вигляді  суми  потенціалів, 

створюваних  окремими  електронами  Фе  =  Хфеї  ,  то  замість  рівняння 

І  =  1 

(28)  одержимо  систему  із  2-2у  рівнянь 


Афеї  =  47іІфіІ^ 


(29) 


де  і  =  1,...,  (2-2у). 
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Перш  ніж  приступити  до  рішення  системи  (29)  необхідно  провести 
усереднювання  |\і/іі^  по  сферичних  координатах  о  і  ф.  Потенціали, 
розраховані  по  приведеній  методиці  для  електронів  від  І8-  до  48-групп, 
представлені  в  табл.  2. 

Таблиця  2 

Потенціали  екрануючих  електронів  остову 


¥ 

Електронні  потенціали 

І8 

Фе  =  е 

Г1  +  2І-І; 

УГ  )  V 

28 

Фе  =е“^'' 

V 

7?  ^  7? 

г  +  г  +  2  + 

8  4  4 

1 

г 

2р 

Фе 

V 

2^  2 

г  +  г  +  2  + 

24  4  4 

1 

г 

38 

Фе  =  е  3 

ґ  4  4  42"^  3  42^  2  42^  5  _  0  1 

е,  2^Г^  +  ^  ^  ^  г+  ^  2+  -  ; 

Тр  З"^  3^  Тр  г 

Зр 

--2г 

Фе  =  Є  3 

Г  2  ^5  4  22^  3  22^  2  42^  5  ^  Л  1 

2^г  +  г-^+  г  +  г+  2+  -  ; 

^3^  2Р  З"^  3^  2р  г  ^  г 

36 

Фе  =  е  ^ 

{  2  ^5  4  ^  2^  3  22^  2  42^  5  ^  0  1 

2^г  +2  +  г  +  г+  2+  -  ; 

81-135  729  81  27  г  г 

48 

-І2г 

Фе  =  е  3 

192^  2 

+  г 

8-192 

^2^  6  52^  5  2  ^5  4  ^  2"^  3 

г'"  -  +  2^г^  -  2  г-^  + 

^4-192^  192^  32-64  8-192 

2^  17^  9^  9 

-  г  +  7+  - 

4  16  4г  ^  4г 

Оскільки  в  моделі  передбачається  розподіл  електронів  по  атомним 
орбіталям,  то  зрозуміло,  що  електрони  дальніх  орбіталей  взаємодіятимуть  з 
меншим  зарядом,  ніж  заряд  ядра.  Тому  при  визначенні  потенціалу  фзі  як 
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величини  заряду  2  необхідно  використовувати  слетерівські  ефективні 
заряди,  тобто  заряди,  одержані  з  урахуванням  константи  екранування  8екр, 
яка  розраховується  за  наступними  правилами: 

•  всі  атомні  орбіталі  класифікуються  по  групах:  (їв);  (28,  2р);  (38, 
Зр);  (36);  (48,  4р)  і  так  далі,  тобто,  8-  та  р-орбіталі  однієї  оболонки 
групуються  разом,  а  6-  та  ґ-орбіталі  представляють  окремі  групи; 

•  константа  8екр  визначається  як  сума  внесків:  0,35  від  кожного 
електрона  в  даній  групі  (для  І8-группьі  0,3);  0,85  від  кожного  електрона 
найближчої  внутрішньої  групи  (в  або  р);  1,0  -  від  кожного  електрона 
більш  глибоких  шарів  (6  або  1); 

•  враховуються  тільки  електрони,  що  знаходяться  між  ядром  і 
досліджуваною  оболонкою.  Електрони,  зовнішні  по  відношенню  до 
досліджуваного,  в  екрануванні  участі  не  приймають; 

•  якщо  досліджуваний  електрон  відноситься  до  рівня  6  або  ґ,  то 
внесок  усіх  внутрішніх  груп  дорівнює  1,0. 

Таким  чином,  8екр  обчислюється  по  наступній  схемі: 

І8  =  0,3-(І8-1); 

28  =  2р  =  0,85-(і8)  +  0,35-(28  +  2р  -  1); 

38  =  Зр  =  1,0-(І8)  +  0,85-(28  +  2р)+0,35-(38  +  Зр  -  1); 

36  =  1,0-(І8  +  28  +  2р  +  38  +  Зр)  +  0,35-(36  -  1); 

48  =  4р  =  1,0-(І8  +  28  +  2р)  +  0,85-(38  +  36  +  Зр)  +  0,35-(48  +  4р  -  1). 

При  обчисленні  констант  екранування  для  даної  електронної 
конфігурації  необхідно  замість  І8,  28,  ...  48  в  правій  частині  рівності 
підставити  число  електронів,  що  заселяють  дану  орбіталь.  Ефективний 
заряд  обчислюється  по  формулі  2еф  =  2  -  8екр.  В  результаті  розрахунків 
одержані  потенціали  фє(г)  електронів  від  І8-  до  48-групп  для  зарядів  2і  = 
1,  22  =  6,  2з  =  28  {рис.  7)  для  деяких  молекул  рудникової  атмосфери. 

Як  приклад  на  рис.  8  показана  залежність  електричного  потенціалу 

від  між'ядерної  відстані  в  молекулі  СО,  створюваного  двома  валентними 
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електронами  в  етані  (1/2,  0,  0)  і  двома  валентними  електронами  в  стані  (3/2, 
0,  0).  Порівнюючи  з  результатами,  представленими  на  рис.  7,  можна 
зробити  висновок,  що  при  розрахунку  електронних  термін  в  першому 
наближенні  екрануючою  дією  валентних  електронів  можна  нехтувати. 
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Рисунок  8  -  Залежність  електричного  потенціалу  від  між  'ядерної 
відстані  в  молекулі  СО:  1  -  валентні  електрони  в  стані  (1/2,  0,  0); 

2  -  валентні  електрони  в  стані  (3/2,  0,  0) 


2.3.  Екранування  ядерних  потенціалів 
валентними  електронами 

Разом  з  екрануванням  потенціалу  ядра  зв'язаними  електронами 
необхідно  врахувати  екранування  валентними  електронами.  Методика, 
покладена  в  оенову  розрахунків,  не  відрізняється  від  наведеної  вище. 

Методику  розрахунку  покажемо  на  прикладі  обчислення 
потенціалу,  створюваного  валентним  електроном  у  стані  (1/2,  0,  0). 
Хвильова  функція  має  вигляд: 

2  2а  •  ехр(-  2а(>^  -  і)) 

11/  =  ^  деа  =  2Ь. 

Х  +  1 
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Розділяючи  змінні,  одержимо 


4к 


дХ 


(х^-1) 


дХ 


^(x)= 


2а  •  е  2а(А.  і) 
Х  +  1 


(ЗО) 


4яА(і_^,2)5ХМ  =  0. 


діл 


діл 


(31) 


Шукатимемо  рішення  рівнянь  методом  варіації  постійних. 
Спочатку  вирішимо  однорідні  рівняння,  внаслідок  чого  матимемо: 


^(x)= 


Сі 


ІП 


>^  +  1 


Х-1 


+  С,, 


(32) 


V(|а) 


2  іл-1 


+  С4, 


(33) 


де  Сі,  Сг,  Сз,  С4  -  постійні  інтегрування. 

Інтегрування  проводимо  за  умови,  що  ф^О  при  нескінченному 
видаленні  від  молекули  і  приймає  кінцеве  значення  при  Х—>1  ц^±1. 
Оскільки  хвильова  функція  фі/2,  о,  о  не  залежить  від  змінної,  то  для 
знаходження  потенціалу  необхідно  вирішити  рівняння  (ЗО). 

Після  нескладних  перетворень  одержуємо 


п 

^2а 


1 

Втг 


Г.  1  ^ 

1п 

Х  +  1 

1-  — 

1  2а; 

1-1 

-2а(^-і) 


+ 


4а^у 


1п 

1-1 


2^Еі(-2а(^  +  і)) 

4а^ 


(34) 

2^Е,(-2а(Х-і)). 

2а  ^ 
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Аналогічно  одержуємо  для  стану  (3/2,  0,  0) 


Фз 


,0,0 


1 

- <: 

Х  +  1 

^-2а(Х-і) 

ґ 

X 

Зтг 

1 

Х-1 

V 

1  + 


13  1 


4а 


( 


а  ^ 


1  + 


1 


1 


7 


.-2а(А,-і) 


2а  4а^  8а^ 

п 

^2а  4а  у 


+  е"^^Еі(-2а(:^-і)) 


^  а  2а^У 
7 


1  З 

—  +  — у  +  — у 

а  4а^  4а^ 


5^1  1 

—  +  — 

а  2а- 


(35) 


Представлені  рішення  задовольняють  умовам  на  нескінченності 
А,^оо  і  при  похідні,  аж  до  другого  порядку,  існують  і  є 

безперервними  функціями. 

Вперше  у  квантовій  механіці  розроблено  метод  розрахунку 
стійкості  речовини,  що  заснований  на  рішенні  задачі  про  рух  електрона 
в  полі  двох  кулонівських  центрів,  що  враховує  вплив  третього 
кулонівського  центру  на  хімічний  зв'язок. 

Рішення  рівняння  Пуасона  (усередненого  по  напрямах  в  просторі), 
знайдені  в  наближенні  центрального  поля  електричного  потенціалу, 
створюваного  зв'язаними  електронами.  Заряди  ядра,  що  входять  в 
одержані  вирази,  розраховувалися  за  правилом  Слетера.  Розроблений 
новий  підхід,  який  враховує  поляризацію  електронних  остовів  при  їх 
взаємному  впливі  один  на  одного.  Разом  з  цим,  одержані  рішення 
рівняння  Пуассона  в  еліпсоїдальній  системі  координат  для  знаходження 
потенціалів,  створюваних  валентними  електронами.  Це  дозволяє 
врахувати  додаткове  екранування  ядерного  потенціалу  електронами,  що 
забезпечують  хімічний  зв'язок  в  молекулі.  Розрахунки  показали,  що  така 
модель  дозволяє  надавати  не  тільки  якісну,  але  й  кількісну  оцінку 
поведінки  дво-  і  триатомних  молекул.  Слід  підкреслити,  що  модель  не 
передбачає  яких-небудь  інших  взаємодій,  а  дає  лише  можливість 

зрозуміти  динаміку  взаємодії  однієї  молекули  з  одним  іоном  певної 
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валентності  [9;  37]. 

Вперше  встановлені  закономірності  поведінки  хімічного  зв'язку  і 
його  енергії  у  присутності  іонів  (де  Е  -  енергія  хімічного 

зв'язку,  К.  -  відстань  між  атомами,  2±  -  знак  і  величина  заряду  іона) 
залежно  від  відстані  між  молекулою  та  іоном,  валентності  іона, 
напруженості  зовнішнього  електричного  поля  та  інших  параметрів  [38- 
39].  Оскільки  реальна  температура  фізичної,  хімічної  та  біологічної 
систем  завжди  перевищує  0  К,  то  не  важко  зрозуміти,  що  із  збільшенням 
температури  функції  Е(2‘^“^;К.)  матимуть  тенденцію  до  більш  сильної 
положистості.  Іншими  словами,  збільшення  температури  системи 
приводить  до  розриву  зв'язків  молекул  на  таких  відстанях  від  точкових 
зарядів,  які  перевищують  у  декілька  разів  характерні  міжатомні  відстані 
у  твердому  тілі. 

Встановлено,  що  відстань  від  іона,  на  якому  відбувається  розрив 
зв'язку,  росте  із  збільшенням  валентності  іона,  густини  розподілу  зарядів 
на  поверхні  або  із  збільшенням  температури  за  інших  рівних  умов  та  із 
зменшенням  енергії  зв'язку  молекули  [9].  Так  у  молекул  СО2  розрив 
зв'язків  відбувається  на  більших  відстанях  від  заряду,  ніж  молекул  СО. 


2.4.  Урахування  впливу  зовнішніх  кулонівських  центрів 
на  електронний  терм  молекули 

При  рішенні  багатьох  задач  квантової  хімії  виникає  проблема 
урахування  впливу  кулонівських  потенціалів  на  електронний  стан 
молекули.  Для  розрахунку  стійкості  кристалічної  структури  речовини 
використовувалася  розглянута  вище  квантово-механічна  модель.  В 
першому  наближенні  поправку  до  електронного  герму  можна  визначити 
з  виразу 
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(40) 


п  N 

АЕ(К)=  ЕЕ 

)=іі=і 


:¥і 


¥і 


^  / 


де  2і/Гц  -  кулонівський  потенціал,  створюваний  і-м  зарядом  2і,  діючим  на 
)-й  електрон  молекули;  н  -  кількість  валентних  електронів,  на  зв'язку. 

Кулонівський  потенціал  в  еліпсоїдній  системі  координат  знову 
представимо  у  вигляді  розкладання  Неймана 
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де  А-<,  А->  -  менша  і  велика  з  координат  X  і  А-і,  Рк'™'  і  ^к'™'  -  функції 
Лежандра  І  і  ІІ-го  роду;  м  =  Х-  = 

Якщо  хвильові  функції  не  залежать  від  координати,  як,  наприклад, 
для  перших  трьох  основних  станів  (15)  -  (17),  то  матричний  елемент 
<\1^^I^VI\1^^>  можна  значно  спростити: 
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Обчисливши  інтеграли  (43),  одержуємо  наступні  вирази  для 
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поправок  до  енергій  електронного  терму  молекули,  що  враховують  дію 
"третіх"  центрів  на  молекулу: 
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На  рис.  9  представлені  розрахунки  електронних  термін  молекули 
N2,  що  знаходиться  в  полі  зовнішнього  кулонівського  центру, 
розташованого  на  різних  відстанях  Н  по  нормалі,  що  проходить  через 
центр  вісі  молекули.  Напруженість  електричного  поля  одновалентного 
іона  на  відстані  10  ^  м  складає  -1,5-10^  В/см.  При  зменшенні  відстані  до 
іону  напруженість  збільшується.  Навіть  при  температурі  0  К 
напруженість  поля  іона  є  руйнівною  для  будь-якого  хімічного  зв’язку. 

Дія  полю  заряду  на  молекулу  призводить  до  "розпушування"  її 
хімічного  зв'язку.  Вплив  зовнішнього  поля  заряду  (кулонівського 
центру)  слабшає  із  збільшенням  відстані  від  молекули.  Закономірності 
змінювання  енергії  зв’язку  молекул  повітряної  атмосфери  О2,  СО2,  та 
СО  від  знаку,  величини  точкового  заряду  і  відстані  від  заряду  мають 
якісну  аналогію  з  поведінкою  зв’язку  молекули  N2  (рис.  9).  Характер 
змінювання  стану  хімічного  зв'язку  мало  відрізняється  від  знаку 
точкового  заряду  (рис.  10). 

2.5.  Визначення  коливальних  рівнів  молекул 

Основним  каналом  розпаду  молекул  і  утворення  атомів  в  плазмі 
кисню  та  азоту  є  дисоціація  прямим  електронним  ударом.  При  аналізі 
механізмів  утворення  атомів  необхідно  враховувати  можливість 
ступеневої  дисоціації  молекул  О2,  СО2  та  N2  при  їх  взаємодії  з 
коливально-збудженими  однойменними  частинками.  Величини  енергій 
розриву  зв’язків  у  молекулах,  як  правило,  значно  перебільшує  величину 
коливального  кванту.  Очевидно,  що  у  процесі  дисоціації  можуть 
прийняти  участь  тільки  молекули,  що  знаходяться  на  високому 
коливально-збудженому  рівні.  У  припущенні  про  триноровський 
розподіл  коливально-збуджених  молекул,  концентрація  збуджених 
частинок  з  енергією,  що  перевищують  поріг  дисоціації,  складає  близько 
10  ...10  см  . 
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Відстань  між  атомами,  К.,  А 


Відстань  між  атомами,  Я,  А 


а 


б 


Рисунок  9  -  Змінювання  потенційної  енергії  Е(К)  молекули  N2  у  полі 
негативного  заряду  X:  а  -  від  величини  2  (відстань  Н  від  заряду 
до  хімічного  зв'язку  2,5-10^^^ м);  б  -  від  відстані  Н  при  2=(^). 

1  -  не  збурена  молекула;  2,  З  та  4  -  відповідно  молекула  у  полі  заряду  2 
(-2),  (-4)  и  (-8);  5,  6  и  7  -молекула  від  заряду  відповідно  на  відстані Н 

(м):  10  \  5-10  ^^  та  1,5-10 


При  знаходженні  коливальних  рівнів  молекули  викориетовуєтьея 
рівняння  Шредінгера,  в  якому  енергія  електронного  терму 
апрокеимується  потенціалом  Морса: 


+ 
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1^0 
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рівноважна  відстань;  м  =  ^1^2  -  приведена  маеа;  А  -  параметр,  що 

визначає  енергію  дисоціації  хімічного  зв'язку;  а  —  знаходиться  з  умови 
перетину  модельної  і  реальної  потенційних  кривих  з  віссю  К. 
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Рисунок  10  -  Характерна  динаміка  змінювання  потенційної  енергії  Е(К) 
хімічного  зв’язку  якої  завгодно  молекули  від  наближення  її  до  точкового 
заряду:  1,4-  молекула  незбурена;  2,  З  -  відповідно  заряд  негативний  (- 
1 )  і  (-2);  5,  6-  відповідно  заряд  позитивний  (+1)  і  (+2) 


Якщо  ввести  нову  змінну  і  =  2л]2МАК^'о  •  !  ах,  то  рівняння  (44) 

перетвориться  до  вигляду: 
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де  введені  наетупні  позначення: 
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Рівняння  (46)  має  два  лінійно  незалежних  рішення 
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Для  утворення  зв'язаного  стану  потрібна  правильна  поведінка 
хвильової  функції  на  нескінченності.  З  цієї  причини  рішенням  може 
бути  тільки  функція  Мп+8+1/2: 
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Якщо  8  і  п  дискретні,  то  ряд  у  фігурних  дужках  обривається.  Це 
приводить  до  умови,  коли  п  =  о,  1,  2...  .  Коливальні  рівні  визначаються 
із  рівняння: 
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Таким  чином,  використовування  (як  наближення  до  електронного 
терму  молекули)  потенціалу  Морса  приводить  до  простого  рішення 
задачі  знаходження  коливального  епектру  молекули. 

На  рис.  11  показані  результати  розрахунку  коливальних  рівнів 
молекули  N0,  що  знаходиться  у  збуджених  станах  (1/2,  О,  0)  і  (3/2,  0,  0). 

При  розрахунку  матричних  елементів  необхідно  провести 
інтеграцію  по  координатах  першого  і  другого  електронів.  Для  цього 
можна  представити  вираз  І/гід  у  вигляді  розкладання  Неймана  в 
еліпсоїдній  системі  координат: 
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велика  або  менша  з  Хі  і  А-г;  Рр'™'(А<)  і  ^р'™'(А>)  приєднані  функції 
Лежандра  І  і  II  роду. 
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де  с  -  постійна  Ейлера. 
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Енергію  електрон-електронної  взаємодії  в  станах  (3/2,  0,  0),  (5/2,  0, 
0)  із  задовільною  точністю  можна  розраховувати  по  формулі  (25),  якщо 
величину  Ь  визначати  з  (14). 

Спроможність  молекули  до  участі  у  хімічних  реакціях 
визначається  рівнем  її  коливального  стану,  тобто  вплив  поля  заряду  на 
хімічний  зв'язок  переводе  молекулу  до  якогось  визначеного  коливально- 
збудженого  стану.  Особливо  це  буде  помітним  під  час  збільшення 
температури  системи:  чим  вища  температура,  тим  на  більшій  відстані 
від  заряду  здійснюється  руйнування  молекули. 


а  б 

Рисунок  11  -  Коливальні  рівні  молекули  N0:  а- в  стані  (1/2,  0,  0); 

б -(3/2,  0,  0) 


Вперше  розроблено  модель  розрахунку  електронних  молекулярних 
термін  з  урахуванням  екранування  ядерних  потенціалів  електронними 
остовами  зв'язаних  електронів.  Апробація  моделі  здійснювалася  при 
оцінюванні  електронних  термін  двоатомних  молекул  СО,  О2,  N2,  ЕіН,  Н2, 
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ЬІ2,  неї  та  інші.  Особлива  увага  приділялася  розробці  методики 
розрахунків  рівнів  коливально-збуджених  станів  молекул  газу  [40].  Ці 
роботи  є  дуже  важливими  при  дослідженнях  хімічних  реакцій,  що 
протікають  у  газах  під  чає  електричного  розряду. 

У  рівнянні  балансу  коливальної  енергії  входять  конетанти 
швидкоеті,  які  описують  релаксацію  молекул  азоту  відповідно  на  N2  і  0. 
Не  дивлячиеь  на  малу  концентрацію  атомів  киеню  у  повітрі  (етаціонарна 
концентрація  киеню  дорівнює  10^^  см~^),  пов’язана  з  ними  релакеація 
носить  неадіабатичний  характер  і  за  швидкіетю  порівняльна  з 
релаксацією  на  азоті.  Швидкості  реакції  та  релакеації  знаходятьея  в 
залежноеті  від  концентрації  атомів  киеню,  які  гинуть  головним  чином  у 
тримолекулярній  реакції  озонування. 

Ветановлено,  що  при  високому  виході  оксидів  суттєвий  вклад 
вносять  також  процеси  доокислення  N0.  Встановлено  також,  що 
утворення  N0  ефективно  здійенюється  у  пасивній  фазі  із  збереженням 
необхідного  ступеню  нерівноважності.  Аналогічні  результати  одержані  і 
для  СО2.  Важливою  загальною  якіеною  особливіетю  нерівноважних 
плазмохімічних  процесів,  що  стимулюються  коливальним  збудженням,  є 
те,  що  чає  реакції  у  пасивній  фазі  набагато  менший  за  час  коливальної 
релаксації. 

Стійкість  оксиду  азоту  в  зворотних  реакціях  визначаєтьея 
головним  чином  за  двома  параметрами:  концентрацією  атомарного 
азоту  і  температурою.  У  першому  випадку  загибель  окиеу  азоту 
незначна  і  може  виконуватиея  лише  за  умов  підвищення  температури 
киеню  за  1000  К.  У  другому  випадку  для  здійенення  зворотних  реакцій 
потрібне  подолання  високого  активаційного  бар’єру  і  лише  при  1500  К. 

З  точки  зору  енергетичної  ефективноеті  найбільш  перспективними 
є  процеси,  що  протікають  селективно  через  коливально-збуджений  стан 
молекул  при  значному  відриві  коливальної  температури  від 
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поступальної.  Такі  процеси,  як  свідчать  теоретичні  розрахунки  та 
результати  експериментальних  досліджень,  найбільш  проето 
реалізуютьея  у  газах  N2,  СО  і  СО2,  які  не  тільки  інтенеивно  коливально- 
збуджуються,  а  ще  й  відносно  повільно  релаксують. 

Змінювання  концентрації  кисню  від  20  до  40%  в  азотно-кисневій 
газовій  суміші  приводе  до  збільшення  розрядної  напруги  на  15%.  Більш 
еуттєві  змінювання  вольт-амперних  характериетик  розряду 
епоетерігаються  при  збільшенні  вміету  озону  в  повітрі. 
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3.  МОДЕЛЬ  КВАНТУВАННЯ  ЕЛЕКТРОНА  В  ПОЛІ 
ЛІНІЙНОГО  ЛАНЦЮЖКА  ІЗ  N  КУЛОНІВСЬКИХ  ЦЕНТРІВ 


Експериментально  встановлено,  що  при  електричному  розряді 
хімічні  реакції  протікають  в  газах  між  негативними,  позитивними  та 
нейтральними  молекулами  і  атомами.  Подолання  енергетичного  бар’єру 
реакцій  полегпіуєтьея  при  учаеті  збуджених  чаетинок.  Хімічні  реакції 
також  протікають  між  поверхневими  атомами  металів  і  збудженими 
молекулами  газу.  В  цьому  випадку  моделювання  імовірних  фізико- 
хімічних  процееів  здійенювалоея  з  урахуванням  оеобливоетей 
ланцюжкових  етруктур. 

Розглянемо  кінцевий  лінійний  ланцюжок  атомів  завдовжки  К. 
Заряд  першого  атома  позначимо  через  2а,  а  заряд  Х-го  через  2ь.  В 
еліпсоїдній  системі  координат  (А-,  ц,  ф) 


.  Га+Гь  . 


в  = 


к 


де  Га  і  Гь  -  відетань  від  електрона  до  першого  і  Х-го  атома  відповідно. 

Розглянемо  математичну  задачу  про  етан  електрона  в  полі  N 
кулонівеьких  центрів.  Для  цієї  моделі  рівняння  Шредінгера  в  атомній 
сиетемі  координат  буде  таким: 


Аф  -  [Е  -  П(2,|іі,ф]-ф  =  0,  (47) 

де  Е  -  енергія  електрона,  ІТ  -  оператор  потенційної  енергії  сиетеми  N 
кулонівеьких  центрів.  Предетавимо  0(2,  ц,  ф)  у  наетупному  виді: 
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■‘а 


+ 


Л,  +  ц  Л,  —  ц 


N-1 

+  ^Фі{7.^,Х,]л,ір)\ 
і=2 


(48) 


Тепер  запишемо  потенціал  Фі  в  декартовій  системі  координат,  при 
цьому  вісь  02  направимо  від  заряду  2а  до  заряду  2ь. 


,  2  2 

ґ 

Кі  ^ 

X  +  у  + 

2  - 

— 

V 

N-0 

(49) 


Далі  робимо  перехід  Фі(2і,  х,  у,  2)  Фі(2і,  X,  ц,  ф),  при  цьому 
використовуємо  наступні  співвідношення 


х^  =  •  (2^  —  і)(і  —  ц^)- со8^  ер;  (50) 

•  (2,2  —  іХі  —  8Іп2  ер;  (51) 


2^  = 


к 


ґ 


Х^  ■  — 


4і 


V 


N-1 


2|іі  + 


4і 


(N  - 1)' 


(52) 


Підставляємо  у  рівняння  (49)  вирази  (50)  -  (52): 


Фі(2і,2,|а,ф) 


22і 

+  |а) 


2^  + 


V 


-1  + 


4і 


N-1 


1 

2 


(К  - 1)2  ^ 
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Після  таких  нескладних  перетворень  потенціал  и(?і,  ц,  ф)  прийме  вид: 


Я 


) - ^ 

А.  +  ц  А.  —  ц 


N-1 


+  Е 

і=2 


2, 


А.  +  ц 


+  Фі(2,ііі,ф,і) 


(53) 


22- 

де  фі  (:а,  р,  ф)  =  ^  •  §і  (:а,  р,  ф,  і), 


§і 


22 


я(а  +  р) 


-1 


ґ 


1- 


(1  +  2Ар)(К  -  іу  +  4іАр(К  - 1)  -  4І 
(А  +  р)^(К-і)^ 


:2  А' 


Тоді  (48)  можна  перепиеати  наступним  чином: 


Щ)і,р,ф)  =  — 


2  2 


+ 


N-1 

+  І 


Я  А-  +  р  А-  -  р 


2  2 

—  +  ц,  Ф,  і) 


А,  +  р  А,  +  р 


(54) 


Знання  хвильової  функції  для  потенціалу  ІІ(А,  р,  ф)  без  фі(А,  р,  ф), 
дає  можливіеть  визначити  хвильову  функцію  для  потенціалу  (54)  по 
теорії  збурення  з  функцією  Гріна  таким  чином: 


Фк,л,п  =  ф(^а .  )  +  1 ї'.  2а ,  2ь ,  2і )  •  со{і',  2а ,  2ь ,  2і )  X 

^  Ф к.Л.п  ’  2а  ?  2^ ,  2^ )  •  бт, 
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ДЄ  Сй(^1 


22; 


+  )л,^ 


2^  +  |л,^  —  1  + 


V 


N-1 


N—1  о  7 


1 


(К  - 1)2  ^ 


N-1 

Якщо  позначимо  2а  =2,  +  2:2;,  2ь  =  2н,  то  одержимо  потенціал, 

і=2 

ЩО  відповідає  рівнянню  (47),  яке  в  еліпсоїдних  координатах  має  вид: 


Р 


0' 


)А,А(і 

'дХ  д[і  ^ 


р 


ар 


+ 


Ф  -  И' 


2: 

їе<р' 


ЛД  + 


(55) 


+ 


2-[е  +  П(?^,р)}д  =  0. 


Для  кулонівського  потенціалу  після  представлення  рішення  у 
вигляді  \\і  =  ^(Л)г|(р)Ф(ф)  рівняння  (55)  розділяється  на  три  звичайні 
диференціальні  рівняння  другого  порядку: 


+ 


Ф(ф)  =  о. 


(1 


Л2 


1-р^ 


—  +  р2  —  К 


л(р)=о. 


А 

дХ 


Л2 

22-1 


+  22є  +  22+  +  К 


де  |ЛІ  -  ціле  число,  є  =  ЕКА;  7г  =  (2а  ±  2ь)-К,  А(К)  -  постійна 
розділення;  Ф(ф)  =  ехр(іЛф). 

Видно,  що  рішення  рівняння  Шредінгера  для  електрона,  що 
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знаходиться  в  полі  N-кулонівськиx  центрів,  розташованих  уздовж 
прямої,  зводиться  до  рішення  вже  відомої  задачі  про  електрон  у  полі 
двох  кулонівеьких  центрів.  Як  приклад  на  рис.  12  приведені  залежноеті 
Е(К.)  електронних  термін  деяких  лінійних  ланцюжкових  молекул. 

В  результаті  фізико-математичного  моделювання  ветановлено,  що 
стійкість  ланцюжка  визначається  його  складом,  співвідношенням  С/Н  та 
енергією  хімічного  зв'язку  С-С.  Міра  етійкоеті  зменшу єтьея  у  такій 
послідовності:  СгНб  ^СзНв  ^С4Ніо  ^С5Ні2,  відповідно  до  якої 
співвідношення  С/Н  збільшуєтьея  (0,333  ^0,375  ^0,400  ^0,416),  а 
енергія  зв'язку  зменшується.  Але  доказів  немає,  що  така  закономірність 
буде  аналогічною  і  для  інших  вуглеводневих  ланцюжків. 


Рисунок  12  -  Електронні  терми  на  прикладі  деяких  вуглеводневих 
ланцюжків:  1  -  С2Н6;  2  -  С3Н8;  З  -  С4Н10;  4  -  СбНд 
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Таким  чином,  надано  теоретичне  обґрунтування  фізико-хімічного 
механізму  впливу  точкових  зарядів  нерівноважної  низькотемпературної 
плазми  і  вторинної  плазми,  утвореної  із  атомів  анода  у  проміжку  між 
електродами,  на  етійкіеть  хімічних  зв’язків  молекул  газу  повітряної 
атмосфери.  Розроблені  теоретичні  основи  динаміки  високопівидкісних 
селективних  хімічних  реакцій,  що  обумовлені  об’ємними  і 
поверхневими  процесами  у  нерівноважній  низькотемпературній  плазмі 
газового  розряду.  Фізико-математичні  моделі  використані  для 
вирішення  проблеми  негативного  впливу  дугогасильних  процесів  в 
оболонках  вибухозахищених  та  рудникових  апаратів,  на  порядок 
підвищити  надійність  роботи  шахтних  та  рудникових  апаратів  захисту 
та  керування  технологічними  комплексами. 

Створена  методика  квантово-механічних  розрахунків  є  новою  у 
квантовій  хімії.  За  допомогою  методики  з'являється  можливість 
одержувати  такі  теоретичні  результати,  що  передбачують  поведінку 
частинок  якої  завгодно  речовини  при  їх  взаємодії.  Тобто  методика  може 
надавати  прогнозну,  але  точну  оцінку  імовірним  фізико-хімічним 
сценаріям,  що  протікають  у  речовині  під  час  тих  чи  інших  збуреннях. 

Заряджені  поверхневі  структури  конденсованих  середовищ, 
кластерні  молекулярні  системи  та  інші  часто  описуються  за 
допомогою  лінійних  атомних  ланцюжків.  У  зв'язку  з  цим  вирішена 
задача  квантування  електрона  в  полі  N-кулонівськиx  центрів.  Знайдені 
перетворення  дозволили  звести  рівняння  Шредінгера  для  ланцюжка  до 
двоцентрової  задачі.  Разом  з  цим,  одержано  аналітичне  рішення  задачі 
про  вплив  третього  кулонівського  центру  на  ланцюжкову  структуру. 

Таким  чином,  основні  фізичні  уявлення  про  елементарні  акти 
хімічних  реакцій  викладені  на  прикладі  реакцій  в  газі.  Проте  подібну 
інформацію  автори  використовували  при  дослідженні  процесів  на 
твердих  поверхнях  [9],  на  межі  розподілу  фаз  [41].  Багато  якісних 
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особливостей  самого  акту  хімічного  перетворення  в  газі  і  в 
конденсованому  середовищі  однакові,  але  в  першому  випадку  вони 
виетупають  в  чистому  вигляді,  не  ускладнені  проблемою  взаємодії 
реакційного  центру  з  середовищем. 

Для  оцінки  етійкості  атомно-молекулярних  систем  розроблено 
фізико-математичну  модель  розрахунку  електронних  термів  молекул  з 
викориетанням  квантово-механічної  моделі  визначення  стану  електрона 
в  полі  двох  кулонівеьких  центрів.  Розроблено  новій  підхід,  що  враховує 
поляризацію  електронних  остовів  при  їх  взаємному  впливі  один  на 
одного,  отримані  рішення  рівняння  Пуассона  в  еліпсоїдальній  системі 
координат  для  находження  потенціалів,  які  утворюються  валентними 
електронами,  що  дозволяє  враховувати  додаткове  екранування  ядерного 
потенціалу  електронами,  що  забезпечують  хімічний  зв’язок  в  молекулі. 

Вперше  були  ветановлені  фізичні  явища,  визначені  їх  якісні  та 
кількісні  параметри,  які  призводили  до  катастрофічно  швидкого 
руйнування  апаратів  управління  й  захисту  на  660  та  1140  В.  Результати 
цих  доеліджень  визначили  шляхи  подальшого  удосконалення 
вибухозахищених  та  рудникових  електричних  апаратів  на  напругу  вище 
500  В.  Були  створені  комутаційні  апарати  з  герметизованими  та  з 
вакуумними  дугогасильними  камерами,  термін  елужби  яких  5... 6  і 
більше  років.  Виробництво  таких  апаратів  в  Україні  було  налагоджено 
на  Торезькому  електротехнічному  заводі  та  підприємстві  „Полярон”  м. 
Львів. 
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висновок 


Встановлено  нове  явище  руйнування  хімічних  зв’язків  молекул  у 
полі  точкового  електричного  заряду,  на  підетаві  якого  запропоновано 
фізичний  механізм  каталітичної  дії  точкового  заряду  та  розроблено 
квантову  динаміку  елементарних  актів,  що  відображають  фізико -хімічні 
особливоеті  виеокопівидкіених  реакцій  в  рудниковій  атмосфері, 
обумовлених  об’ємними  та  поверхневими  процееами  у  нерівноважній 
низькотемпературній  плазмі  газового  розряду. 

Надано  теоретичне  обґрунтування  фізичного  механізму  впливу 
точкових  зарядів  нерівноважної  низькотемпературної  і  вторинної 
плазми  у  проміжку  між  електродами  на  етійкість  хімічних  зв’язків 
молекул,  що  дозволило  вирішити  проблему  негативного  впливу 
дугогаеильних  процееів  в  вибухозахищених  оболонках  і  на  порядок 
підвищити  надійність  роботи  шахтних  та  рудникових  апаратів  захисту 
та  управління  технологічними  комплекеами. 

Вперше  в  облаеті  квантової  хімії  розроблено  фізико-математичну 
модель  поведінки  хімічного  зв'язку  молекули  в  полі  точкового 
електричного  заряду  -  кулонівського  центру  та  модель  динаміки 
елементарного  акту  хімічної  реакції  у  газі  під  час  впливу  зовнішнього 
поля.  Основну  теоретичну  трудність  представляло  рішення  задачі  про 
поведінку  двоатомної  молекули  в  полі  кулонівського  центру  (іону 
металу)  тому  у  даній  роботі  були  використані  результати  наших 
попередніх  досліджень  в  облаеті  кінетики  утворення  нової  фази.  Деякі 
основні  особливості  вирішення  задач  запропонованої  моделі  наетупні. 

Рішення  рівняння  Пуасона  (усередненого  по  напрямах  в  просторі), 
знайдені  в  наближенні  центрального  поля  електричного  потенціалу, 
створюваного  зв'язаними  електронами.  Заряди  ядра,  що  входять  в 
одержані  вирази,  розраховувалиея  за  правилом  Слетера.  Розроблений 
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новин  ПІДХІД,  який  враховує  поляризацію  електронних  остовів  при  їх 
взаємному  впливі  один  на  одного.  Разом  з  цим,  одержані  рішення 
рівняння  Пуассона  в  еліпсоїдальній  системі  координат  для  знаходження 
потенціалів,  створюваних  валентними  електронами.  Це  дозволяє 
врахувати  додаткове  екранування  ядерного  потенціалу  електронами,  що 
забезпечують  хімічний  зв'язок  в  молекулі.  Розрахунки  показали,  що  така 
модель  дозволяє  надавати  не  тільки  якісну,  але  й  кількісну  оцінку 
поведінки  дво-  і  триатомних  молекул.  Слід  підкреслити,  що  модель  не 
передбачає  яких-небудь  інших  взаємодій,  а  дає  лише  можливість 
зрозуміти  динаміку  елементарного  акту  хімічної  реакції,  тобто  взаємодії 
однієї  молекули  з  одним  іоном  певної  валентності. 

Вперше  встановлені  закономірності  поведінки  хімічного  зв'язку  і 
його  енергії  у  присутності  іонів  (де  Е  -  енергія  хімічного 

зв'язку,  К.  -  відстань  між  атомами,  2±  -  знак  і  величина  заряду  іона) 
залежно  від  відстані  між  молекулою  та  іоном,  валентності  іона, 
напруженості  зовнішнього  електричного  поля  та  інших  параметрів. 
Оскільки  реальна  температура  фізичної,  хімічної  та  біологічної  систем 
завжди  перевищує  0  К,  то  не  важко  зрозуміти,  що  із  збільшенням 
температури  функції  Е(2‘^“^;К.)  матимуть  тенденцію  до  більш  сильної 
положистості.  Тобто  збільшення  температури  системи  приводить  до 
розриву  зв'язків  молекул  на  таких  відстанях  від  точкових  зарядів,  які 
перевищують  у  рази  характерні  міжатомні  відстані  у  твердому  тілі. 

Встановлено,  що  відстань  від  іона,  на  якому  відбувається  розрив 
зв'язку,  росте  із  збільшенням  валентності  іона,  густини  розподілу  зарядів 
на  поверхні  або  із  збільшенням  температури  за  інших  рівних  умов  та  із 
зменшенням  енергії  зв'язку  молекули.  Так  у  молекул  СО2  розрив  зв'язків 
відбувається  на  більших  відстанях  від  заряду  ніж  молекул  СО. 

Розроблена  модель  розрахунку  електронних  молекулярних  термін  з 
урахуванням  екранування  ядерних  потенціалів  електронними  остовами 
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зв'язаних  електронів.  Апробація  моделі  здійснювалася  при  оцінюванні 
електронних  термін  двоатомних  молекул  СО,  О2,  N2,  ЬіН,  Н2,  ЬІ2,  НС1  та 
інші.  Особлива  увага  приділялася  розробці  методики  розрахунків  рівнів 
коливально-збуджених  станів  молекул  газу.  Ці  роботи  є  дуже 
важливими  при  дослідженнях  хімічних  реакцій,  що  протікають  у  повітрі 
під  час  електричного  розряду  в  газах. 

Відкрита  нова  закономірність  розвитку  високошвидкісних  фізико- 
хімічних  реакцій  в  рудниковій  атмосфері,  обумовлена  об'ємними  і 
поверхневими  процесами  у  нерівноважній  низькотемпературній  плазмі 
газового  розряду,  що  супроводжуються  руйнуванням  молекул  та 
утворенням  нових  сполук  під  впливом  електричного  поля  точкових 
зарядів  як  каталізаторів  хімічних  процесів. 

Оскільки  осередки  плазми  під  час  розряду  виникають  не  тільки  у 
газі,  а  і  на  металевих  поверхнях  катоду  вирішувалося  питання  відносно 
сценарію  поверхневих  хімічних  реакцій.  Встановлено,  що  для  розриву 
хімічного  зв'язку  необхідною  і  достатньою  умовою  є  наявність 
точкового  заряду,  дія  його  поля  на  молекулу  в  процесі  зближення 
молекули  із  зарядом  на  відстань,  при  якій  зв'язок  молекули  руйнується. 
В  цілому  ж  активність  хімічних  реакцій  на  поверхнях  росте  із 
збільшенням  температури  системи  за  рахунок  збільшення  амплітуди 
коливань  поверхневих  атомів. 

Розроблений  квантово-механічний  метод  розрахунку  стійкості 
речовини,  заснований  на  рішенні  задачі  про  рух  електрона  в  полі  двох 
кулонівських  центрів,  що  враховує  вплив  третього  кулонівського  центру 
на  хімічний  зв'язок.  Метод  апробований  при  розрахунку  стійкості 
лінійних  ланцюжків  атомів  водню  і  вуглецю.  Одержані  електронні  терми 
і  міжатомні  відстані  задовільно  узгоджуються  з  експериментальними 
даними.  Рішення  квантово-механічної  двоцентрової  задачі  узагальнено 
на  N-кулонівськиx  центри.  Цей  підхід  дозволив  одержати  рішення  задачі 
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квантування  електрона  в  полі  лінійного  молекулярного  ланцюжка. 
Зокрема,  показано,  що  при  збільшенні  кількості  валентних  електронів 
підвищується  мінімум  енергії  ланцюжка  і  збільшуються  міжцентрові 
відстані.  При  певній  кількості  «зайвих»  електронів  ланцюжок  стає 
нестійким. 

За  даними  офіційних  джерел  інформації  з  кожним  роком  зростає 
негативний  вплив  результатів  діяльності  людини  на  стан  навколишнього 
середовища,  заглиблюються  встановлені  раніше  проблеми  та 
ускладнюється  їх  вирішення,  збільшується  число  і  концентрація  викидів 
антропогенних  речовин  у  атмосферу  і  таке  інше.  Значна  частина 
атмосферних  викидів,  що  обумовлена  горінням  палива,  надходить  від 
теплових  електростанцій.  Внесок  у  забруднення  атмосфери  із  щорічним 
збільшенням  темпів  роблять  підприємства  хімічної,  металургійної  та 
нафтопереробної  галузей.  Головними  компонентами  згоряння  та 
переробки  хімічної  продукції  є  оксиди  вуглецю,  сірки,  азоту. 
Запропонована  модель  може  бути  корисною  у  проектуванні  запобіжних 
заходів  щодо  зниження  навантаження  на  екологічний  стан 
навколишнього  середовища  антропогенних  речовин. 
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